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Controlling Fringe Sensitivity of Electro-Optic 
Holography Systems Using Laser Diode Current 

Modulation 

ABSTRACT 

Electro-Optic Holography (EOH) is a non-intrusive, laser-based, displacement 
measurement technique capable ol static and dynamic displacement measuiements. EOH 
is an optical interference technique in which fringe patterns that represent displacement 
contour maps are generated. At excessively large displacements the hinge density may 
be so great that individual fringes are not resolvable using typical EOH techniques. This 
thesis focuses on the development and implementation of a method for controlling the 
sensitivity of the EOH system. This method is known as Frequency Translated Electro- 
Optic Holography (FTEOH). It was determined that by modulating the current source of 
the laser diode at integer multiples of the object vibration, the fringe pattern is governed 
by higher order Bessel function of the first kind and the number of fringes that represent a 
given displacement can be controlled. The reduction of fringes is theoretically unlimited 
but physically limited by the frequency bandwidth of the signal generator, providing 
modulation to the laser diode. Although this research technique has been verified 
theoretically and experimentally in this thesis, due to the current laser diode capabilities it 
is a tedious and time consuming process to acquire data using the FTEOH technique. 
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CHAPTER 1.0 


INTRODUCTION 

Holographic interferometry has been used as a method for determining suitace 
displacements for small amplitude vibrations for a number of years. This laser-based 
technique is extremely sensitive and is capable of measuring displacements as small as 

— where X is the laser wavelength. The interferometric method uses the principles ot 
100 

light interference to create a fringe pattern to measure structuial displacement [1]. Using 
this method for out-of-plane measurements one complete (black to white) interference 
fringe is approximately equal to one-half laser wavelength of light [2]. The fringe 
patterns are captured and analyzed producing displacement information at a sensitivity 

level that would otherwise be unattainable. 

The goals undertaken by the author of this work are laid out in Section 1 .3 and 
include the method for controlling EOH fringe density and conversion of the current 
mirror based system to a fiber optic based system. In Section 1.4 and 1.5 background 
information on the physics and behavior of light as well as an introduction to holography 
is presented. 
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1 1 NOMENCLATURE 


CCD Charged Couple Device camera 

BOH Electro-Optic Holography 

FTEOH Frequency Translated Electro-Optic Holography 

k Wave number, equal to 2 ni X 

X Wavelength of the laser light used to produce interference 

6 The angle between the illumination and reflection beams 

1.2 PROBLEM STATEMENT 

Conventional interferometry, although sensitive and relatively accurate, has many 
difficulties which hinder the rate of data acquisition. Conventional holography requires a 
photographic recording and development process preventing real-time measuiements. It 
is this limitation which instigated the development of a more advanced method for 
performing holographic interferometry. Using this method interference patterns are 
formed and recorded with a charged coupled device (CCD) camera and processed by a 
computer acquisition system. The technique is known as Electro-Optic Holography 
(EOH). This method is used to generate fringe patterns in a near real time fashion with 
the same sensitivity as typical holographic interferometry [3], 

The major drawback of EOH is that the measurement range is limited. For 
example, a typical laser diode with a wavelength of 852 nm has an upper measurement 
limit of approximately 20 microns. With increased displacement, the number of fringes 
that appear in the interference pattern increase causing the EOH image to become 
saturated with fringes. When the fringe density becomes too high it is impossible to 
clearly identify individual fringes, which prevents counting fringes accuiately. This lack 
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of resolution is the physical mechanism limiting the upper measurement range of the 
EOH system. In order for large amplitude displacements to he measured the sensitivity 
of the EOH system must be reduced such that each fringe represents a displacement 

, A 

greater than — . 

C.C. Aleksoff [5] discussed a conventional holographic interferometry technique, 
which involved modulation of the reference beam as means ot controlling fringe 
characteristics. He stated that by 
modulating the reference wave at 
higher harmonics any number of 
fringes corresponding to the smaller 
amplitude displacements could be 
eliminated, thereby alleviating the 
image of high fringe density. This 
method, however, has yet to be 
established with EOH technology. 

This thesis describes the method for 
controlling fringe density for EOH measurements, thereby allowing larger displacement 
measurements to be acquired. The goal is to take an image like that of Fig. 1 . 1 and reduce 
the number of fringes so that resulting the image will look like that of Fig 1 .2. 

A simple flow chart showing the modem methods for determining structural 
displacements is shown in Figure 1.3. The figure illustrates where FTEOH will be placed 
within the discipline of current metrology. The flow chart is broken down into two 
general categories: theoretical/numerical methods and experimental methods for 



Figure 1.1: An EOH 
image for an object 
under dynamic loading, 
too many fringes for 
clear resolution. 



Figure 1.2: An 
optimal EOH 
image, with 
resolvable fringes. 
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determining structural displacements. hOH can be described as a subset of holography 
under the experimental methods category. HOH itself can be broken down into two 
additional sub-categories. Typical HOH and HTHOH, which is the conti ibution this thesis 
work will add to the arena of structural displacement measurement methods. 

In Chapter 2 the EOH theory is presented in detail and mathematical models for 
the intensity and phase distribution for both the static and dynamic cases are examined. 
Finally, Chapter 3 includes the method for controlling fringe density for the EOH system 
by means of laser diode current modulation and the results are verified. 
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Figure 1.3 Thesis contributions to the field of EOH. 
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1 3 THESIS GOALS 


• Develop a method for controlling EOH fringe density using 
laser diode current modulation. 

• Convert the existing optical setup from a mirror-based setup to 
an optical fiber based setup. 

• Modification of the computer software to allow for FTEOH data 
analysis. 
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1 4 INTRODUCTION TO LIGHT 


Light has proven to be of major interest by many great scholars. Newton, 
Huygens, Young, and Rayleigh, rather than shown light could behave dually, attempted 
to prove that light was either a particle or a wave [6], In the seventeenth century Isaac 
Newton was the most prominent advocate ol the particle theory. Newton had proclaimed 
that rays of light were streams of particles, which were emitted bom a souice and 
traveled in a straight line. His supporting argument stated that light from a source casts a 
sharp shadow of objects in its path, in contrast to water and sound waves, which bend 
around obstacles [7], therefore light must be a particle. 

In contradiction to Newton and also a contemporary of the same time was Dutch 
physicist Christian Huygens who favored the opposite theory that light was a wave. 
Huygens’ discussed the wave behavior of light by explaining that light spreads out from a 
light source in all directions not just in a straight path. Thomas Young, an Englishman in 
support of Huygens, performed an experiment in which he passed light through two small 
slits and produced a complex interference pattern of light and dark bands a phenomenon 
that could only be explained by the fact that light must be a wave. 

The beginning of the nineteenth century served to show strong support for the 
wave theory but as the century came to a close so too did this support [7], Such physicists 
as Planck, Einstein, Bohr, and Compton once again proved that light must be a particle by 
examining blackbody radiation, the photoelectric effect and the scattering of X-rays from 
electrons [7], 

In 1924, Luis de Broglie began to look at light with a new approach when he 
began to look at the dichotomy of light. He suggested that a particle with momentum had 
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associated wave properties. This new idea brought the argument ol light full circle but 
now under the pretense of duality. James Clerk Maxwell was able to combine these 
concepts and determined that light is merely electromagnetic energy, which he described 
as electromagnetic waves [8]. Maxwell determined that the complete spectrum of 
electromagnetic energy includes radio, infrared, visible, ultiaviolet, x-ray and gamma 
radiation. His theory also encompassed the quantum theory, as revealed by Planck, 
Einstein, Bohr, and Compton, who stated that electromagnetic energy is ‘quantized,’ and 
can be taken from the electromagnetic field in discrete amounts known as photons. Thus, 

the modern concept of light was bom [8], 

Although light does have a dual nature, the wave behavior holds the greatest 
interest in the study and research of Electro-Optic Holography, which is the focus of this 
thesis. Therefore, a general understanding of wave properties as applied to EOH will be 
examined in detail. The inquiry will be inclusive to Young’s double slit experiment, 
diffraction, wave interference, and superposition, which are the foundations of the EOH 
technique. 

Thomas Young provided the first demonstration of the wave characteristic of 
light in 1801, by demonstrating that light exhibits an interference behavior. He 
established that light waves coming from a single source arriving at a point by two 
different paths can combine and cancel each other by destructive interference to produce 
an interference pattern as shown in Fig. 1.4. Constructive interference occurs when the 
slopes of the two interfering waves carrying the same sign add together to create a band 
of light in the interference pattern. Conversely, destructive interference occurs when the 


8 


slopes of the two waves are opposite in sign and when added togethei cancel one another 
out to create a band of dark in the interference pattern. 


Wave 



Figure 1.4: Construction and Destruction of Light waves as described by Young. 

At the time the interference phenomenon was discovered, it was not explainable 
by means of the particle theory. It was thought that two or more particles could not come 
together and cancel one another out in such an orderly pattern. Therefore, the 
phenomenon had to be the result of wave interaction. 

Consequently, the wave theory was not immediately accepted among scholars. 
There were those who understood that waves, such as sound, traveled through some type 
of medium. Their argument was based on the fact that light from the sun arrives at the 
surface of the earth through a vacuum of empty space, without a medium. Therefore, 
light must be a particle because waves need a medium tor propagation [6]. In the mid- 
nineteenth century, British physicist James Clerk Maxwell suggested that a medium, 
which he called the ether, must exist everywhere, even in the empty vacuum ot space and 
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it was this substance, he proclaimed, that light, as a wave was able to propagate. 
However, in 1887, Albert Michelson and Edward Morley at the Case school ot Applied 
Science in Cleveland, Ohio, disproved the existence of ether. Dutch physicist Hendrtc 
Lawerence and Albert Einstein [9] later disproved the existence ol ether thereby leaving 
the wave theory in limbo. 

Nevertheless, the argument proceeded on with one of the most intriguing 


questions of the time: If the nature of light we 
to bend around objects like water and 
sound waves thus allowing one to see 
around comers. Of course we know that we 
cannot actually see around corners, but in 
fact light does bend around them. Foi 
example, if an object is held up against a 
light so that a shadow is cast, one should 
notice that the edges of the shadow are not 
as sharp as one might expect. These soft 
edges are a result of light bending around 
the object. This type of bending is known 
as diffraction. Diffraction is the modulation 
of waves in response to an obstacle or 
grating in the path of propagation. Figure 
1.5 is an illustration of Young’s double slit 
experiment, which demonstrates diffraction 


e to be wave-like, then light would be able 
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Figure 1.5: Young’s Double Slit Experiment 
|6], (A) Light traveling through a large 
diameter hole does not experience diffraction. 
(B) Light traveling through a small diameter 
hole with respect to the wavelength is 
diffracted upon reaching exit of the hole. 
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on slit openings of different sizes. Diffraction of light will occur regardless of the size of 
object, but in Fig. 1.5 (A) the effect of diffraction is less apparent because the size of the 
slit opening is relatively large with respect to the wavelength of the illuminating light. 
Therefore, the light rays that meet the wall of Fig. 1.5 (A) will continue to move through 
the opening in a straight line. On the other hand, if the diameter of the hole is small with 
respect to the wavelength of light, like it is shown in Fig. 1.5 (B), the effect of diffraction 
is more apparent. The light passing through the relatively small aperture will spread out 
in all directions as a result of diffraction. 

When coherent light waves pass through a small aperture like that ol Fig. 1.5 (B), 
an interference pattern, rather than a sharp spot of light is seen on the second screen [6], 
Because there are two slits in Young’s experiment, the light from the two slits interfere 
with each other in both a constructive and destructive manner, which is dependent on 
path length. Therefore, diffraction can be regarded as a consequence of interference from 
many coherent wave sources. In other words, diffraction and interference are basically 
equivalent [6]. 

The short wavelengths of visible light make interference patterns difficult to see, 
and in order to observe this phenomenon the following conditions must exist: 

1 The light sources must be coherent and of identical wavelength. 

2. The superposition principle must apply. 

The first condition states that the two interfering waves must be coherent, by 
definition, longitudinal coherence is a function of laser frequency stability and emission 
linewidth, not phase. However, one can argue that two incoherent beams cannot maintain 
a constant phase relationship since they will be of different frequencies. Two waves are 
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said to be in phase with one another when the wave crests and valleys ot both waves 
coincide. If all three waves of Fig. 1.6 coincide, then they are in phase with one anothei 
regardless of amplitude. However, two waves are said to be out ol phase when the wave 
crests and valleys do not coincide but instead intersect elsewhere along the wave. Thus 
the frequency along the wave must be maintained, otherwise, the wave is said to be 

incoherent [6]. 



Figure 1.6: Light waves in phase . 

The second condition states that the superposition principle must apply, which is 
simply the arithmetic addition of the two interfering waves. Two coherent waves 
originating from the same source only differ by the path length for which the two waves 
travel. The difference in path length causes the two waves to differ in optical phase. This 

difference, upon superposition, results in interference. 

As mentioned previously, Thomas Young was the first to demonstrate the wave 
behavior of light through means of interference in his double slit experiment, as shown in 
Fig. 1.5. A drawing of the complete experiment is shown in Fig. 1.7 (A) with the 
resulting interference pattern Fig. 1.7 (B). In the experiment, the light is incident on a 
screen in which there is a small opening [6], The light emerges from the opening in the 
screen and is incident on a second screen that contains two small openings. These two 
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openings serve as a pair of coherent light sources because the waves emerging from them 
are created by the same wavefront and, therefore, maintain a constant phase relationship. 
When the light reaches the third screen the interference pattern ol Fig. 1.7 (B) can be seen 
[6], The areas of constructive interference are described in Fig. 1.7 (A) as ‘max and 
result in a band of light in Fig. 1.7 (B). The opposite it true for destructive interference, 
which will produce both a bright fringe (constructive interference) and a dark fringe 
(destructive interference). This phenomenon is illustrated in Fig. 1.8. The two waves 
leaving the two small openings leave in phase and strike the third screen at a central point 
in Fig. 1.8 (A). At this 
point the two waves have 
traveled an equal distance. 

The waves are said to be in 
phase, thus resulting in 
constructive interference or 
a band of light in the 

interference pattern. In Fig. 

, , Figure 1.7: (A) Young’s Double Slit experiment and 

1 .8 (B) the path lengths are the ^ interference pattern produced |6|. 

not equal because the top wave has to travel one wavelength further than the bottom 
wave. Since the path difference is one wavelength the two waves are still in phase and 
will again result in a band of light. Last, in Fig. 1.8 (C) the top wave has traveled one half 
wavelength farther than the bottom wave. Therefore, the wave valley of the top wave will 
add to the wave crest of the bottom wave causing destructive interference or a dark band 
in the optical interference pattern. 
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(A) 


( B) 


(C) 


Figure 1.8: (A) The light from the two holes travels equal distance and reach the screen in phase. 
(B) The light from the top hole travels one wavelength further than the light lrom the bottom 
hole, light is in phase. (C) The light from the top hole travels one half of a wavelength further 
than the liuht from the bottom hole, the licht is out of phase |6|. 


The light and dark bands in the interference pattern in Fig. 1.7 (B) appear to be 
speckled with light and dark spots throughout. This phenomenon is known as speckle 
noise and is a result of light scattering off of an optically rough surface. The roughness of 
the surface varies the optical path length causing constructive and destructive interference 
with the diffracted wave front, which results in a granular, speckle appearance thioughout 


the interference pattern. 


1.5 INTRODUCTION TO HOLOGRAPHY 


Many optical measuring techniques take advantage of light interference using the 
fact that light and dark bands, fringes, represent approximately one wavelength of light. 
This technique is known as interferometry. One such interferometric method is known as 
holographic interferometry. The images produced from this method are known as 
holograms, which will now be discussed. 

In 1948, British scientist Dennis Gabor invented a photographic technique known 
as holography. Holography produces a three-dimensional image that has all the depth and 
parallax of a real-life scene [7], Gabor called his images holograms and explained their 
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three-dimensionality as a reconstruction of the coherent scattered light leflected off of a 
structure combined with a coherent reference light field [7-10]. Optical coherence, or the 
ability of a light source to maintain constant phase, is ol vital importance in the 
development of holograms. In 1948, however, a coherent light source was not readily 
available, and the ability to produce holograms was at best difficult. The growth of 
holography stood at a virtual stand still until 1962, when coherent lasers became 
available. The advent of the laser brought back interest to holography, and Emmet Leith 
and Juris Upatnieks of the University of Michigan produced the first laser based 
holographic image [7], Leith and Upatnieks' development resulted in a surge of activity 
with scientists world wide trying their hand at holography and the production of 
holograms [11]. 

A hologram is a three dimensional image ol an object produced by recording the 
patterns of interference, on a special photographic film, formed by a split laser beam. 
The holographic image can only be viewed after "reconstruction. Reconstruction is the 
illumination the photographic plate with the same wavelength lasei light, illuminated at 
the same "reference beam angle". In the case of a simple hologram the image is created 
by first splitting the beam of a coherent laser into two separate beams, the object beam 
and the reference beam . The object beam is used to illuminate the test object, and the 
second beam is maintained as a reference source. Upon illumination of the test object, 
the beam is phase modulated, with respect to the reference beam [12]. The modulated 
beam is recombined with the reference beam to yield an interference pattern. The 
interference pattern can then be seen on the special photographic film when illuminated 
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by the reference beam. figure 1.9 illustrates a simple optical set up used foi creating a 


holographic image. 



Figure 1.9: The basic set up for creating a hologram. 

The resultant interference pattern ox fringe pattern includes the phase 
relationships of each part of the light from the scene with the reference light and therefore 
with every other part [7], The interferogram is essentially a coded contour map where the 
interference fringes are the lines of code that contain the information about the depth and 
displacement of the object. The displacement information can be decoded by counting the 
number of fringes which represent lines of constant displacement, multiplying the count 


by the quantity 


X 


, where 6 is the angle between the illumination and reflection 

2k cos 6 

vectors and k is the wave number (Section 2.1 details these vectors and derives the above 
quantity in full). The method of determining displacement through the use of 


interference fringes is called holographic interferometry. 
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Holographic interferometry provides the means by which deformations and 
displacements can be analyzed through the use of a interferogram. Holograms have 
traditionally been recorded on light sensitive photographic 1dm that operates much like 
diffraction gratings [13]. During interferometric applications a large number of 
holograms are needed in order to fully analyze the displacements of the object of interest. 
After the holographic images have been recorded on film, the interpretation and 
quantitative analysis of the interferograms begins. This process can be tedious and 
extremely lengthy, especially when dealing with interferograms of dynamically excited 
objects. More importantly this recording and development process inhibits real-time 
measurements. 

Since the development of holographic interferometry a more advanced method for 
processing holographic images has arisen. This method records interfeience patterns 
with a charged coupled device (CCD) camera and analyzes the data with a coniputei 
acquisition system. The three-dimensionality of the images is lost using this new 
technique but the rate of computation is much faster. This technique is EOH and can 
generate interference patterns in a near real time fashion with the same sensitivity as 
typical holographic interferometry [3]. 
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CHAPTER 2.0 


EI.ECTRO-OPTIC HOLOGRAPHY 

Holograms have traditionally been recorded on either photographic him or 
expensive cumbersome glass film [13]. When large numbers of holograms are required, 
inexpensive, flexible film is usually chosen. Regardless, film is difficult to work with. 
The development of EOH began in 1985, presenting a method to capture and stoie 
holograms electronically [13]. It was during this time that Karl Stetson [13] and others 
[14-16] came to develop EOH. This chapter discusses EOH in detail with a brief 
introduction in section 2.2 followed by a detailed description of both static EOH and 
dynamic EOH in sections 2.3 and 2.4 respectively. The chapter will conclude with a 
discussion of the sensitivity vector. 

2.1 NOMENCLATURE 

A 0 Amplitude of the object beam, the laser beam directed onto the test object. 

A r Amplitude of the reference beam, the laser beam directed into the CCD 

camera. 

B Phasor bias, a phase change incurred by the reference beam to solve for £2, 
the fringe locus function. 

CCD Charged Couple Device camera 

e„ The amplitude of vibration for the test structure 
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E„ The time dependent displacement o( a point on the vibrating stiucture. 

EOi 1 Electro-Optic Holography. 

c, The amplitude of the PZT. 

E, The time dependent displacement of a point on the vibrating structure. 

/' Emission frequency of the laser diode. 

F interference The resulting light field of the scattered object and reference beams. 

F o The object light field, the beam of light that illuminates the test object. 

F r The reference light field, the beam of light that is directed into the CCD 

camera. 

FTEOH Frequency Translated Electro-Optic Holography. 

F 1CC The digitally reconstructed light field. 

F v The resulting object beam upon reflecting off of the test object. 

< 7 (J A geometrical factor with = 2 for displacements perpendicular to the 
object surface with normal incidence. 

1 The light intensity of the laser beam. 

l al The time averaged background intensity. 

I mt The time dependent maximum intensity. 

l n The intensity of the n th sequential frame. 

T n The intensity of the n th sequential frame with the object in its deformed state. 

I m The time dependent intensity of the n’ h sequential frame. 

I () The intensity of the object beam, the beam directed onto the test object. 

F 0 The intensity of the object beam while the object is in its deformed state. 

I,„ The time dependent intensity of the object beam. 

l r The intensity of the reference beam, the beam directed into the CCD camera. 
J„ The zero order Bessel function of the first kind. 
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Sensitivity Vector, the vectorial difference between the propagation vectois. 
The wave number whose value is given by 2 k/ A,. 

Vector of propagation for both illumination and leflection, whtic / — 1 to 4. 
Magnitude of displacement of the test object. 

Distance the light travels. 

The characteristic equation given by the zero order Bessel function of the 
first kind. 

The fringe order. 

The point location of the light source. 

The point of interest on the test object. 

The point of interest on the test object while the object is in its displaced 
state. 

Piezoelectric actuator. 

Point detector. 

Direction vector from the point source O to the point detector Q. 

Direction vectors where / = 1 to 4. 

Exposure time of the CCD camera. 

Displacement in the z-direction. 

The optical phase difference between the scattered object beam and the 
reference beam. The magnitude is an unknown value. 

The time dependent unknown phase difference between the scattered object 
beam and the reference beam. 

90° sequential phase step used to determine the fringe locus function. 

Error between the estimated fringe locus function and the actual fringe locus 
function. 



0 Arbitrary phase value assigned to the object and scattered object beams. 

0 Optical phase of the object beam. 

0 Optical phase of the reference beam. 

0 Optical phase of the illuminating object beam. 

0 ^ Optical phase of the reflected object beam. 

0 Angle between the illumination and reflection vectors. 

A The wavelength of the laser light used for interference. 

£!(/) The fringe locus function, 

n The approximated fringe locus function. 

2.2 Introduction to EOH 

EOH is a non-intrusive, laser-based, whole-field measurement technique that uses 
the basic physics of holographic interferometry to measure small surface displacements 
on statically or dynamically loaded structures. The fringe pattern that results from the 
optical interference of scattered object light combined with reference light are essentially 
contour maps of the object’s displacement, where each fringe represents a line of constant 
displacement [17]. EOH differs from typical holography in that the images are not 
recorded on special photographic film and processed individually but instead are captured 
with a charged-coupled device (CCD) camera and processed with a computer acquisition 
system. The captured images are digitized and processed all within a matter of seconds to 
produce a near real time measurement [18]. This method allows the deformation of the 
object’s surface to be directly perceived by the observation of the fringes [19], where a 
great number of fringes would indicate a relatively large displacement. 
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Phase 

Stepper 


Laser Wavelength = 852 nm 

Figure 2.1: Optical set-up |18|. 


HOH interference patterns are generated by the combination of scattered light 
from the test object and non-scattered light from a reference source. This is 
accomplished by splitting the beam of a coherent, single frequency laser into object and 
reference beams. The object beam maintains most of the laser light energy, while the 
reference beam is typically of about 5% of the laser energy [4], As shown in Fig. 2.1 a 
microscope objective is used to expand the object beam and project it onto the test 
specimen. A CCD camera lens focused on the test object collects a portion of the 
scattered light and directs it to a beam combiner, where it is recombined with the 
reference light. The collected scattered light has a random, unknown phase distribution, 
which cannot be measured directly by the CCD camera. However, 
constructive/destructive interference caused by relative phase differences between the 
collected scattered object light and reference light will result in a spatial intensity 
variation dependent upon the relative phase. The CCD camera records these spatial 
intensity variations in the form of an electronic hologram. Movement of the test object 
results in changes of the object beam optical path length, producing a change in the 
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collected scattered light/reference beam relative phase and obseived pixel intensity [4]. 
For the case of static measurements FOH interference patterns, or interjerograms , result 
when two electronic images of the object obtained under ditferent states of stress aie 
digitally subtracted. In this manner, differences in optical path length between the non- 
loaded and loaded conditions can be distinguished as intensity variations observed by the 
CCD camera [4]. The method of subtracting the non-stressed image with the stressed 
image is known as double exposure holography and is used to determine displacements 
for statically loaded structures. 

2.2.1 EOH Theory. The Dynamic Case 

For the case of dynamically loaded structures, which is the majoi focus of this 
thesis, the double exposure method is not sufficient in gathering all the information 
necessary to obtain displacement information [20], To do so, a method known as 
continuous time averaged interferometry is used. This method requires that an 
interferogram of the moving object be recorded for an exposure time well in excess of the 
period of oscillation. By guaranteeing that the exposure time is greater than the period of 
oscillation, the interferogram will effectively record a continuous distribution of double 
exposure EOH images corresponding to the various points during the vibration cycle 
[20]. This interferogram represents the time-averaged displacement of the object as well 
as the intensity variations of the interference fringes, for the period of oscillation. Just as 
the displacement varies with the vibration cycle so does the intensity distribution of the 
interference fringes. It can be shown that the intensity of the interference fringes varies as 
the square of the zero order Bessel function of the first kind. This relation can be derived 
by first defining the light fields from both the object F,„ and reference F„ beams as well 
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as the light field scattered from the vibrating object I\ Next, by defining the interaction 
between these light fields, the Bessel function as it relates to the intensity distribution of a 
dynamically loaded object, can be derived as follows. 

[■' ( v, v./ ) = A it cxp[i’0 H (x,y./) + iii(t) + 2;z// ] (2.1 ) 

where A„ is the amplitude of the object beam, lift) is the temporal phase change caused 
by reflection off the vibrating structure, / is the emission frequency of the laser diode and 
is the phase of the object beam given by Eq. 2.2 




2xg u E a (x,y,t) 

A 


( 2 . 2 ) 


where g„ is a geometrical factor with g„ = 2 for displacements perpendicular to the object 
surface with normal incidence. E„ is the time dependent displacement of a point on the 
vibrating structure and can be expressed as 

E i ,(x,y,t) = e ll sm(cot). ( 2 - 3 ) 

In Eq. 2.3 e„ is the amplitude of vibration for the test structure and (O is the frequency of 
vibration. The reference beam field can be expressed as Eq. 2.4 
F r (x,y,t) = A r exp[/0„(x,y,O + 2nf't] 


where A, is the amplitude of the reference beam and <p r is the phase, which is given as Eq. 


2.5 


Q,Ax,y,t) 


2^g,.£ r (x,>-,t) 

A 


(2.5) 


where g r is a geometrical factor with g, = 2 for displacements perpendicular to the object 
surface with nonnal incidence. E, is the time dependent displacement of a point on the 
vibrating structure and can be expressed as 
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v,0 = e r sin(ax) . 


(2.6) 


In Eq. 2.6 c, is the amplitude of the PZT and (0 is the frequency of vibration. 

When the reflected object beam recombines with the reference beam, the phase 
difference between the two beams form a pattern of fringes. These fringes are described 
by the fringe locus function £2(t), which is the major element loi determining stiuctural 
displacements. However, in order to solve for the fringe locus function, the intensity 
relation of the fringes must first be solved. Interference is the addition of the refeience 
and the reflected object beams, which can be expressed as Eq. 2.7 

F , = [A exp (i(4>„ + £2(0 + 2 reft) + A, exp (<j), + 2 tft )] (2.7) 

The intensity, or the energy delivered by the wave per unit area per unit time [6], is 
defined as the product of the interference beam (Eq. 2.7) and its complex conjugate, 
which can be rewritten as Eq. 2.8 

/~[F r + FJF r + Fj (18) 

where * represents a complex conjugate. Substitution of Eq 2.1 and 2.4 into Eq. 2.8 
yields the intensity relation of the interference fringes as shown in Eq. 2. 9- 2.22. 

/ = [A o exp (i<f> a + /£2 + 2 rft)+ A,, exp (i<j>, + 2 rft )]• ^.9) 

[A 0 exp(- i{<t> 0 + £2 + 2 jtft))+ A r exp- i{<p r + Irft )] 

Recognizing that exp(i'ry) = cos <y + sin ft) and exp(-/w) = cosft)-sin<y and 

expi{27tft -27rft) = \ , Eq. 2.9 can be expanded as shown in Eq. 2.10 

/ = [4 (cos((Z> ( , + £2) + i'sin(0 () + £2))+ A r (cos^ ; + isin^>, )]• 

[/4 ( ,(cos^„ +Q)-/sin(0„ +£2)+ A r cos<p, — / sin 0, ] (2.10) 

Multiplying the dot product of Eq. 2.10 and using the identity cos .v + sin x — 1 and 

i 2 = —1 (complex variable) Eq. 2.10 expands to the following: 
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/ = /!,; + A‘ 2 + A„A r [ cos(0„ + S2)cos0, - /sin <p t cos(0„ + Si)] 

+ 4 4 ]/sin(0„ -t-S2)cos^ + sin($, + £2)sin(9 i ] 

+ A„A\ cos<£, cos( 0„ + 12)-isiii(0„ + £2)cos$ ] 

+ 4 , 4 , [/sin0, cos(<Z> + S2) + sm0, sin(d>, + 12 )} ( 2 . 1 1 ) 

Eq. 2.1 1 can be reduced as shown by Eq. 2.12 

/ = 4,;4, 2 + 24„4 ( [cos(0„ + 12)cos0, +sin(d>„ + S2)sin0, ]. (2.12) 

Recognizing the trigonometric identity cos( A' — y) = cos .v cos y + sin a sin y . Eq. 2.12 is 
reduced to the following 

/ = A* A, 2 + 24, ,4,. cos[(0„ + i1)-<p, ]• (2.1o) 


Where F V F’ = F r F‘ and is given by 

F F' = [4„(cos(0„ + Q) + tsin(0,))][4, (cos(0, ) — !sin(^>, ))]. 
Multiplication of Eq. 2.14 yields 
= 4„4 r cos((0„ +Q.) -(/>,.), 

hence. 


(2.14) 


(2.15) 


(2.16) 


FF; + F r F‘ = 2cos((& + 12) -(/>,)■ 

The intensity relation can now be rewritten as Eq. 2. 1 7 

/ = A] + A 2 r +FF; + F r F' . (217) 

The camera will only record the average energy over the duration of the exposuie time. 
Thus, Eq. 2.17 must be integrated over time in order to obtain the correct intensity 
relation. This operation yields Eq. 2.18 

1 „ r' . 1 C 1 j. (2.18) 


T 


J Idt = 4 + A;+- F r £ Kdt + - f; J (1 Fdt . 
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The final term in Eq. 2.14 is of primary importance; it governs the form of the EOH 
interference fringes. Thus by evaluating this term it can be deteimmcd that the intensity 
of the fringes will be directly proportional to the zero order Bessel function as shown in 
Eq. 2.19-2.22. 

— f ' Fdt = — fVexpfiW,, +il(t)))dt < 219 ) 

'[' JO 1 J Jo ' 

When the EOH image is digitally reconstructed, it follows that the reconstructed 
wavefront F rcr is given by Eq. 2.20. 


F = — A cxpf/f ) f exp(/Q(/ ))<:// (2.20) 

tv< y ° Jl) 

The integral of Eq. 2.20 is called the characteristic function and is denoted by M as 
shown in Eq. 2.21 


^£exp(/Q(0)^ = M(i}(0) (221) 

The time average integral of Eq. 2.21 can be evaluated as T 0 (Q) . The zero order Bessel 
function of the first kind with argument O(f) and the intensity of the reconstructed 
wavefront will then be proportional to Eq. 2.22 


/“LC = /( » ex P(' 2 ^) 


— I exp(/Q)£fr 




( 2 . 22 ) 


where, q is a constant. Thus, from the physics of two interfering beams off the surface of 
a vibrating object the intensity of the fringe is found to be proportional to the zero order 
Bessel function as shown in Eq. 2.22. 
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2.2.2 Verification of the Bessel Function Intensity Distribution 


The zero order Bessel function and its square are shown in Fig. 2.3. As shown in 
Eq. 2.22, the EOH fringe spatial intensity distribution is proportional to the zero order 
Bessel function for vibrating structures. Nodal positions on the stiuetuie can be 
identified by what is termed the zero order fringe, which corresponds to the largest peak 
of J-„. The intensity distribution of this fringe is approximately 6 times brighter than all 
other fringes, providing easy visual identification. The zero order fringe is the brightest 
fringe because of the time integration properties of the camera and scattering statistics 
[4], The zero-order fringe occurs at the nodal positions, where there is no object motion. 
Since there is no object motion at the nodes, pixel intensity is constant and does not vary 
with time. Integrating this over the exposure time of the CCD camera is simply 
integrating a constant and the result is the constant itself. Thus, when two holograms of a 
vibrating object are subtracted, the non-nodal pixels will not be as bright as the zero- 
order fringe [4], When there is an optical phase difference the intensity of the fringes 
decreases until the two light fields are 180° out of phase. This event corresponds to the 
Bessel function’s touching of the x-axis, which results in a dark fringe within the EOH 
interferogram. 
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Figure 2.2: EOH 
image of a 
dynamically loaded 
structure. The large 
bright fringes are 
denoted as zero 
order fringes. 


J<i and Jo 


lo 


lb 


Figure 2.3: Zero order and square of the zero order Bessel function 


Figures 2.4 and Fig. 2.5 illustrate a cantilever beam vibrating at different 
frequencies. In Fig. 2.4 the beam is vibrating such that two zero order fringes appear 
corresponding to the first and second flexural modes. Figure 2.5 illustrates a third 
flexural mode shape for this cantilever beam. These figures are shown as an example of 
how EOH images corresponds to object displacement and motion. Note how the intensity 
varies across the image, where high intensity relates to small displacements and low 
intensity relates to large displacements. 
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Figure 2.4: Second mode shape and corresponding hologram |3| 



Figure 2.5: Third mode shape and corresponding hologram |3| 
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2.3 FRINGE LOCUS AND THE SENSITIVITY VECTOR 


Electro-optic holographic numerical analysis, in particular the detei urination ot 
the fringe locus function £2(0, depends on knowledge of the illumination and the 
observation directions used while recording the EOH images. These directions are given 
by the illumination vector k t and the observation vector A , , which describe the 
propagation of light from a point source to the structure, and from that point, to the CCD 
camera. The difference of these two vectors is known as the sensitivity vector and is 

denoted as K = (k 2 —k i ), where A. is considered to be a vector conesponding to the 

measurement direction which bisects the 
illumination and observation vectors of a 
single point. The fringe locus function is 
related to the sensitivity vector as well as 
the displacement vector L , as shown in Eq. 

2.23. 

KL = Q. (2.23) 

The derivation of the sensitivity vector can 

easily be determined by examining the 

static measurement case, however, the 

vector is the same for both static and 
dynamic cases. First assume the test object is a simple cantilever beam moving in the 
transverse direction as shown in Fig. 2.6 (a). A wave traveling normal to the plane of the 



Figure 2.6: Normal deformation of a cantilever 
beam: (A) Side view of the beam with 
displacement in the z direction and illumination 
in the -z direction, arrows indicate direction of 
illumination. (B) Interference pattern observed 
on the front of the beam |13|. 
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object illuminates its surface, as shown in Fig. 2.6, where the arrows noimal to the 
surface indicate the direction of illumination [22], For the static case, the light scattered 
by the test object is initially recorded while the object is in a no-load condition. Then a 
second capture is identically made with the object under loaded conditions. Interference 
patterns, like the one in Fig. 2.6 (B), are produced when the difteience is taken between 
the light scattered from the no-load conditions and the loaded conditions. The figure 
shows what the virtual image of the reconstruction may look like. The ioice applied to 
the test object is small and displacement only occurs in the z direction, which is denoted 
as z = Z(x). If the distance the light travels is /„ then the distance it would travel after 

displacement is / 0 -2 Z(jc) and the corresponding optical phase shift is equal to 

Q = r^2Z(jt) [22]. The next step in determining displacements is to assign fringe order 

A 

to the bright fringes. If the base of the test object is fixed, the zero order fringe will 
appear at the base of the structure and will be assigned the number N = 0. Sequential 
bright fringes are then assigned positive integers from N =1, 2, 3, . . . N. The A/ bright 

fringe will correspond to a phase change of Q = 2nN = y 2Z(x) . Thus solving for the 
displacement Z, one would arrive at Eq. 2.24 



Displacements normal to the surface of the object are easily determined by counting the 
fringes. However, for structures that undergo more complex motion, (i.e., two and three 
dimensions) the method of determining displacement by simply counting the fringes is no 
longer valid. Vector displacements must be extracted from the interference pattern by a 
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careful point-by-point analysis. First, a model must be developed from which the 
derivation of fringe interpretation can be based. Consider the surtace ol the test object as 
a collection of point scatterers. When a point P is displaced by a distance ol T to its new 
position P' the interference between the light scattered at P and the light scatteied at P 
contribute to the fringes in the interferogram. Figure 2.7 shows the important elements of 
the physical system [22]. 




Figure 2.7: Nomenclature for fringe analysis. (A) Schematic diagram of the system. (B) Position 
and propagation vectors [22|. 




In Fig. 2.7 (B) several vectors are used to determine the relationship between the 
fringe locus function and the displacement vector. The vectors R and r he m the same 
plane as the points O. P and Q. The vectors A, and A, are the propagation vectors ol the 
light illuminating point P from point O and the light rellected point P to point Q 
respectively [22], 

The phase of the illumination beam and the reflected scatteied light beam would 
reach the detector at point Q and would have values as shown by Eqs. 2.25 and 2.26, 

0 =k r r l +k 1 -(R-r l ) + j u (225) 


and 


0, = A, • r, + A 4 (/? -r,) + 0 (( 


(2.26) 


where 0, is given as the arbitrary phase assigned to the beams from the point O and 0 is 
the phase of the scattered light before displacement. The variable 0 is the phase ol the 
scattered light after displacement. The total phase of the scattered beam would then 
reach point Q with a value equal to Eq 2.27 

<2 ' 27) 

After the test object is displaced, the illumination vector and reflected scattered light 
vector are A, and A, respectively. Combining these vectors into Eq. 2.27 the fringe 

locus function can now be expressed as Eq. 2.28 

Q = (k 2 -A,)(/j — F 3 ) + (A, — A,)-r, + - k 2 ) (R - r 2 ) . (2-28) 

In a real EOH system the magnitude of the vectors r t and r, are much larger the 

displacement vector L = |r, - ?j|, so for all practical purposes (A, - A, ) is perpendicular to 
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r, and (A, - A : 2 ) is perpendicular to (R -F,). This condition causes the dot product of 
these four terms to be identically equal to zero thus allowing the fringe locus function to 


(2.29) 


become 

Q = (A 2 -*,)•! 

The difference between these two vectors is known as the sensitivity vector, which is 
denoted by K = (As -A,). The fringe locus function can now be rewritten as 

o-JM. (230) 

In Fig. 2.7 (A) the angle between the propagation vectors is given as 26 [_2]. The 
magnitude of the sensitivity vector can be determined using vector algebra and the unit 
directional components of the propagation vectors as shown by hqs. 2.31 and 2.32 

A : , = |A, | sin(0)t + |A, |cos(0) / ^ 2 ‘ 31 ^ 

and 

As = |A 2 |sin(0)»-|* 2 |cos(0)y. (232) 

The sensitivity vector then becomes the difference of these two vectors given as 

K = 2a(a 2 - A,)cos(0) (2 33) 

where k = 2n / A is the wave number. A single point observation of the test object would 

yield a displacement shown by Eq. 2.34 

r _ Q (2.34) 

2k cos(0) 

The ideal EOH system is one where the propagation vectors lie normal to the plane of the 
test object's surface, and the angle between the vectors is 6 = 0 . However, alignment of 
the equipment usually does not allow for such ideal conditions, and the vectors do have 
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some small angle. The angle is said to be small enough that it is approximately equal to 
zero, which allows the displacement vectoi to become 

J = £_ (2.35) 

2k ' 

where cos(O) = 1 . The full three-dimensional displacement vector is given as 

I = L i + L J + L z k . < 2 ‘ 36) 

CCD cameras used in the EOH technique have a two-dimensional array of 
sensors, or pixels. Each pixel on the CCD array functions as an independent 
interferometer. The light within the speckle field has a particular amplitude and phase 
relative to the reference beam [19]. This amplitude and phase are encoded as intensity 
variations detected by the CCD. By recording and processing the values of the light 
intensity before and after the test object moves, an estimate ot surface deformation can be 
calculated by subtraction of the phase measured before and after loading of the object 
[19]. Although each point will have a different intensity that vanes sinusoidally 
according to its initial phase, the phase change at different points will be the same as long 
as the surface displacement is the same. The result of subtraction of the images at 
different times is the contour map representing areas of constant displacement by means 
of interference fringes [19]. 

Referring back to Eqs. 2.25 and 2.26, the phase variations of a single point are 
based on the fact that when coherent light illuminates the optically rough surface of a test 
object, variations in the roughness of the surface produce diffracted waves. These 
diffracted waves can interfere with each other when the optical path lengths between 
surface variations differ. These local fluctuations cause the image to have a granular 
appearance known as speckle. However, these variations in intensity are directly related 
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to variations in phase, which also vary accordingly across the tough surface [22], If the 
alignment of the EOH system is correct, then the scattered object and reletence beams 
will be coherent with one another and any small variation in optical path length will 
cause the speckle pattern to “sparkle." This indicates that the speckle pattern is moving 
in and out of phase, which allows one to know when the alignment of the EOH system is 
correct. If this phenomenon does not occur prior to data collection, fringes will not be 
formed because the conditions for interference will not have been met. 

2.4 STATIC MEASUREMENTS 

EOH can be broken down into two categories, static measuiements and dynamic 
measurements. Static measurements are obtained by utilizing the double exposure 
method on stationary objects subject to static loading. Dynamic measurements employ 
the time-average holographic method and are used to obtain vibration amplitude data. 
The method of measuring a dynamically loaded object is of major interest tor controlling 
fringe characteristics, the subject of this thesis. In order to completely understand the 
physics of the dynamic case, the technique for obtaining static measurements must first 
be understood. 

2,4,1 STATIC EOH THEORY 

Static measurements are made by capturing two sets of four sequential frames of 
EOH images. Each image varies from the other by introducing a 90° optical phase step 
induced on the reference beam. After the original laser beam has passed through the 
beam splitter, the reference beam phase is directed onto a piezoelectrically driven mirror. 
The PZT moves in such a manner that the reference beam phase is incremented by 90 
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for each frame [3]. The first four images are ol the object in its original position while the 
second four images are of the object in the displaced position. The following equations in 
this section and in Section 2.3 are used for the static and dynamic measurements 
respectively. The equations are used lor the quantitative deteimination of the 
displacement vector and are a product of the work shown by R. Pryputniewicz in 
Holottraphic Numerical Analysis [3], 

The image of the object in its original, unstressed (no displacement) state can be 
described by the illumination of the object in its original state of the nth sequential frame, 
as shown in Eq. 2.37. 

l„(x,y) = lo + Ir - 2 Ao A, cos[A<p + A0„]. (2 -37) 

The illumination or irradiance of the object in the displaced state can be desciibed as 

l,/(x,y) = 1,,' + I.- + 2 A,,' A, cos[Atp + Q - A0„]. (2.38) 

In the above equations, I„ is the intensity of the light scattered off of the object prior to 
deformation, I,,' is the intensity of the scattered light off of the object after deformation 
and I, is the intensity of the reference beam. A0 n is the 90° sequential phase step, which is 
induced onto the reference beam between recordings of each individual frame. The 
addition of the phase step A0 n is used to eliminate the phase difference A<p between the 
object and reference beams, which is an unknown, from the intensity. These steps are 
used to solve for Q, the fringe locus function, which is directly related to the 
displacement of the object. Therefore, adding in the sequential phase steps into Eqs. 2.37 
and 2.38 the four sequential frames recorded by the CCD camera for both positions of the 
object become 

1 1 = lo + l r + 2 A„ A r cos[A<p + 0°], (2.39) 
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1 : = I„ + I r + 2 A„ A, cos[A(p + 90°], (2.40) 

h = I„ + l r + 2 A„ A, cos[A<p - 1 80°]. (241) 

U= |„ + 1,-2 A„ A, cos[Atp - 270°]. (2-42) 

and 

1/ = |,/ + I,. + 2 A,,' A, cos[A(p + £2 + 0°], (2-43) 

I { = 1„' + I r + 2 A,,' A, cos[Acp + £2 + 90°], (2-44) 

1,' = I,/ + i r + 2 A,/ A, cos[Atp - £2 + 1 80°], (2-45) 

U = I,/ + Ir + 2 A,/ A r cos[A<p + £2 + 270°]. (2-46) 


Using both sets of intensity equations, prior to deformation (Eqs. 2.39-2.42) and after 
deformation (Eqs. 2.43-2.46). there are four unknowns - the irradiance I 0 , 1,/, the phase 
difference between the object light field and the reference light field Acp, and the fringe 
locus function, Q. By adding in the 90° phase step A0„. Aq> can be eliminated from the 


intensity equations as shown in Eqs. 2.47 -2.55. 

1 , = I„ + I r + 2 A 0 A r cos(A<p + 0°), (2-47) 

1 2 = I„ + l r + 2 Ao Ar sin(Atp + 90°), ( 2 -48) 

1 3 = 1 () + I r - 2 A 0 A r cos(Acp +180°), (2-49) 

u = I« + Ir - 2 Ao A r sin(Atp + 270°), (2 50) 

and 

I,' = I,/ + 1/ + 2 A 0 ' A r cos(A<p + £2 + 0°), (2-51) 

l 2 ' = i o ' + 1/ + 2 A„' A r sin(A<p + E2 + 90°), (2-52) 

1/ = f/ + 1/ - 2 A 0 ' A r cos(Atp + £2 + 1 80°), (2-53) 

1/ = 1,,' + 1/ - 2 Ao' A r sin(Acp + Q + 270°). (2-54) 
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Using the trigonometric addition and subtraction formula oiT.q. 2.55 


cos(x + y) = cos x cos y sin x sin y (2.55) 

Eqs. 2.47-2.54 become 

1 1 = 1„ + I r + 2 A„ A, cos(Acp), (2.56) 

>2 = I«> + I, + 2 A„ A, sin(Atp), (2.57) 

1, = 1„ + I r -2 A 0 A r cos(Acp), t 2 - 5 *) 

I 4 = I () + I r - 2 A„ A r sin(Acp), (2 59) 

and 

I,' = 1,/ + 1/ + 2 An" A r cos(Atp + Q ), (2.60) 

W = 1 0 ' + 1/ + 2 A,/ A r sin(A<p + £1), (2-6 1 ) 

W = i; + 1/-2 A,,' A, cos(Acp + Q), (2-62) 

u' == I 0 ' + 1/ - 2 A„ A r sin(A(p + £2). (2.63) 


Subtracting Eqs. 2.56 and 2.58 as well as Eqs. 2.57 and 2.59 the following set of 
equations are obtained for the intensity difference image of the object in its undeformed 
state. 

(Ii -1 3 ) = 4 A 0 A r cos Acp (2.64) 

and 

(I 2 - U) = 4 A 0 A r sin Acp (2-65) 

Performing the same procedure for Eqs. 2.60 and 2.61 as well as 2.61 and 2.63 the 
intensity difference images of the object in its deformed state are obtained as the 
following 

(h' - Ij') = 4 A,/ A r cos A<p (2.66) 

and 
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(I 2 ' - I/) = 4 A,,' A, sin Atp 


(2.67) 


The addition of Eqs. 2.64 and 2.65 will yield Eq. 2.68 

(1, 1 3 ) + (1,' 1 3 ') = 4 A,, A, cos Atp * A,/ A, cos(Atp + £2) (2-68) 

Using the trigonometric addition formula cos (A<p T £2) = cos Atp cos£2 - sinAtp sin£2 and 
assuming that the amplitude of the original object beam is approximately equal to the 
amplitude of the displaced object beam, A„ = A,,', and Eq. 2.68 can be rewritten as 

(|, i 3 ) + (I,' - | 3 ') = 4 A„ A, [cos A(p + cos Atp cos£2 - sinAtp sin£2] (2.69) 

and then simplified as Eq. 2.70 

(|, i 3 ) + (1,' - I,') = 4 A„ A, [(1+ cos £2) cosAtp - sin Acp sin £2], (2.70) 

Equation 2.68 is the entire data image processed by the HOH computer. The remaining 
images are then processed in order to solve for £2, the fringe locus function. 

Similarly for Eqs. 2.66 and 2.67 the image of intensity can be given as Eq. 2.7 1 

(I, - 1 4 ) + (l 2 '- V) = 4A 0 A r [(l+ cos £2) sinAtp - cosAtp sin £2], (2.71) 

Finally, Atp is eliminated by adding the squares of Eqs. 2.70 and 2.71 together 

[(li - E) + (Ii' - b')] 2 + [(b - I 4 ) + - UO] = 

{4 A 0 A r [(l+ cos £2)cosAtp-sinAtp sin £2]} 2 

+ [4 A 0 A,[(1+ cos £2) sinAtp-cosAtp sin £2]} 2 (2.72) 

Let [(I, - l 3 ) + (I,' - b')] 2 + [(b - I 4 ) + (b' - b')] 2 be denoted as D, 2 + D 2 2 , then solve to 

eliminate Atp as shown in the following equations 

D| 2 + D 2 2 = 16 A u 2 A, 2 [(1+ cos £2) cosAtp - sinAtp sin £2] 2 

+ 16 A 0 2 A r 2 [(1+ cos £2) sinA-cosAtp sin £2] (2.73) 

Expansion of Eq. 2.73 yields Eq. 2.74 


41 



(2.74) 


D, : + d 2 2 = 16 A 0 ‘ A, 2 [(cos 2 Aip + siirAcp) + 2 cos £2(cos"Aip + sin"A<p) 
+ (cos 2 A(p cos 2 £2 * sm : A<p sin" £2) 

4 - (cos 2 A(p sill 2 £2 + sm 2 A(p cos" <.2)]. 


Recognizing trigonometric identities Eq. 2.74 can be rewritten as Eq. 2.75 

D, 2 + D 2 2 = 16 A, 2 A r 2 1 1 4- 2 cos£2 + [cos 2 (Acp - 12) + sin 2 (Acp - 12)]! ■ (2.75) 

Using the relation (sin 2 x + cos 2 x) = 1 Eq. 2.75 reduces to the form shown in Eq. 2.76 

D| 2 + D 2 2 = 16 A,, 2 A, 2 (1-2 cos£2 r 1). (2.76) 

Rewriting Eq. 2.76 the solution for Dr + D 2 " becomes Bq. 2.77 

D| 2 + D 2 2 = 32 A, 2 A, 2 ( 1 + cosi2) . < 2 - 77 


Using the half angle identity of Eq. 2.78 

(1 + cos £2) = 2 cos 2 (Q/2), ( 2 - 78 ) 

Eq. 2.77 can be reduced to Eq. 2.79 

[D| 2 + D 2 2 ] l/2 = 8 A 0 A, cos(E2/2). (2.79) 

Equation 2.79 represents the static viewing image. In order to obtain the displacement 
data for the EOH system the fringe locus function Q must be determined. In order to do 
so the same procedure used to obtain Eq. 2.77 must be followed for the subtraction of 
Eqs. 2.64 and 2.67 as well as Eqs. 2.65 and 2.67 as shown in Eqs. 2.80 and 2.81. 

(Ii - 1 3 ) - (h' - V) = 4 A 0 A r [( 1 - cos Q.) cosAtp + sinAtp sin Q] (2.80) 


(I 2 — Lj) - (I/ - I 3 ') = 4 A () A r [( 1 - cos £2) sinAip - cosAip sin £2] (2.81) 

Equations 2.80 and 2.81 are then squared and added together as shown in Eq. 2.82 

[(I, - I3) - (1/ - I/)] 2 + [(h - 14 ) - (V - h')f = 

1 6A 0 2 A r 2 [( 1 -cos£2)cosA<|) + sinAcp sinQ] 2 
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(2.82) 


- 16A 0 ‘A r 2 [( 1 - cos S 2 ) sinAtp - eosAtpsin S2 ] 2 . 


Let [(I, 


h)-(b' L')] : ~[(L U ) - ( 1 1 ' 


1 ,')] : be denoted as Dr - Dr and Hq. 2.82 can 


then be expanded to hq 2.83 


Dr + Dr = 1 6 A ,, 2 A r 2 {[cos 2 Atp + sin : Atp] 2 cos' Atp cosi2 

+ [cos 2 Atp cos!2 + siirAtp sinS2] - 2 sin'Atp cosS2 
+ siirAtp eos'S 2 + cos'Atp sin'S2J (2.83) 

Recognizing the identity (cos 2 Atp + sin 2 Atp) = 1 and (cos 2 AtpcosQ 4 - siirAtp sinQ) = 

cos 2 (Atp - S2) Eq. 2.83 reduces to Eq. 2.84 

Dr + D 4 2 = 16A„ 2 A, 2 { 1 - 2 cos 2 Atp cosS2 - cos'(Atp - S2) - 2 sin Atp cosS2 

+ sin 2 A(p cos : S2 + cos 2 Atp sin : S2 \ . (2.84) 

Grouping similar terms Eq. 2.84 can be rewritten as Eq. 2.85 

D 3 2 + D 4 2 = 16A 0 2 A, 2 { 1 + (-2 cos 2 Atp cos£2 - 2 sin' Atp cosE2) 

+ cos 2 (Atp - E2) + (cos 2 Atp sin 2 S2 + sin 2 Atp cos 2 Q){. (2.85) 
Using the identity (cos 2 Atp sin 2 Q + sin 2 Atp cos 2 S2) = sin 2 (Atp -Q) Eq. 2.85 reduces to Eq. 


D 3 2 + D 4 2 = 16A„ 2 A r 2 { 1- 2 cosQ (cos 2 Atp + sin 2 Atp) 

+ cos 2 (Atp - £2) + sin 2 (Atp - £>)} . ( 2 - 86 ) 

Next, using (sin 2 x + cos 2 x) = 1 Eq. 2.86 becomes 

D , 2 + D 4 2 = 1 6A () 2 A r 2 {1-2 cos£2 + 1 } , (2 ‘ 87) 


which reduces to Eq. 2.88 


43 



D 3 2 + D 4 2 = 32A„ 2 A r 2 (l - cosQ). 


( 2 . 88 ) 


Equation 2.77 is subtracted from 2.88 to yield Eq. 2.89 

D = (D | 2 + D 2 2 ) + (D 3 2 + D 4 2 ) = 64 A,, 2 A , 2 cos£2 (2.89) 

Thus far there are two unknowns £2 and D and only one equation, so an additional 
equation is needed in order to solve ior £2. The second equation will be formed by 
following the same procedure as was used to determine D but instead ot working with the 
combinations of Eqs. 2.64 with 2.66 and Eqs. 2.65 with 2.67 the combination will be 
switched to Eqs. 2.64 with 2.67 and Eqs. 2.66 with 2.65. This procedure is as follows 


(li - b) + (I2' - b') = 4 A o 2 A ' 2 [( 1 + sin ^)cos A( P + sinA( P cosQ l ( 2 - 9 °) 

(E - 1 4 ) + (1 1 ' - I./) = 4 A u 2 A , 2 [(1 + sin Q)sin Atp - cosAtp cos£2] (2.91) 

and 

(I, _ i 3 ) . ( 1 / _ l 4 ') = 4 A ,, 2 A , 2 [( 1 - sin Q)cos Atp - sinAtp cos£2], (2.92) 

(1 2 - U) - (1/ - h') = 4 A „ 2 A r 2 [( 1 - sin £2)sm A<p - cosAtp cos£2], (2.93) 

Where [ (I, - I 3 ) + (I2' - V )] 2 will be denoted as N | 2 and [(b - I4) + (b' ~ b')]' Wl11 be 
denoted as N 2 2 . The sum of these two terms yields 

N | 2 + N 2 2 = 32 A 0 2 A r 2 (1 + sin £2). ( 2 - 94 ) 

Let [([, _ i 3 ) . ( 1 / _ L ,')] 2 be denoted as N 3 2 and [(1 2 - I 4 ) - (1/ - b ')] 2 be denoted as N 4 2 
and the sum of these two terms yields 

N 3 2 + N 4 2 = 32 A 0 2 A r 2 (1 - sin £2) (2.95) 

These two terms are then subtracted and denoted by the letter N for simplicity 

N = (N , 2 + N 2 2 ) + (N 3 2 + N 4 2 ) = 64 A ,, 2 A r 2 sin £2. (2.96) 

Dividing Eqs. 2.96 by 2.97 yields Eq. 2.97 

N 64A 2 A 2 sin£2 (2.97) 

D ~ 64A;Afcos£2 ‘ 
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Recognizing the trigonometric identity tun Q — sin 12 / cos Q, Eq. 2.97 can be solved for 
Q, the fringe locus function in terms of N and D as shown in hq. 2.98 

Q = tan' l (N/D). ( 2 - 98 ) 

Finally, the displacement data can be determined using Eq. 2.30 as shown in Section 2.3. 

2.5 DYNAMIC MEASUREMENTS 


This section provides a detailed mathematical derivation of the physical and 
analytical processes required to obtain quantitative vibrational amplitude measurements 
using EOH. Dynamic measurements are made using the time-averaged hologram method 
discussed in section 2.2, where Eq. 2.13 gives the fringe intensity distiibution. The 
intensity can be rewritten as Eq. 2.37 and multiplied by the characteristic equation of Eq. 
2.21 to yield Eq. 2.99. 

I nt (x,y) = lot + Irt + 2 Act A r cos[A(p, + A9 n ]M(£2 t ). (2.99) 

The images of intensity are gathered in a series of sequential frames as shown by Eq. 
2.99, where the n indicates the n h sequential frame [3], which varies as a function of time 
given as the subscript t. M is the characteristic function and is the modulating term of Eq. 
2.99, which is evaluated as the zero-order Bessel function of the first kind. 

2.5.1 DYNAMIC EOH THEORY 

Just as in the static case, Eq. 2.99 contains four unknowns I Ul , , A<p, and £2, . 

In order to determine the fringe locus function from the interferogram of a vibrating 
object, four sequential frames are recorded with the phase ot the reference beam 
incremented in steps of 90° between each frame. This process is represented by the 
following set of simultaneous equations: 
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/„ = /,„ + /,. + 2 A ol A r cos( &<p + 0 " )M (£2, ) 

(2.100) 

l y = / i+ 1+2.4 it A r cos(A0 + 90" )M (12, ) 

(2.101) 

/ ( , = 1'" + /,. + 2A„A, cos( A<p + 1 HO" )M(L1, ) 

(2.102) 

/ M = I ni + /, + 2A M A r cos(A0 -t- 270" }M(Q, ) . 

(2.103) 

In the above equations the argument (x, y) has been dropped tor 

simplification and the 

subscript t indicates the time varying parameters. Following the same procedure used tor 

static measurements the equations are reduced to 


/„ = /„, + + 2 A ol A, cos(A$ )M (Q, ) , 

(2.104) 

/ : , = + A + 2 4„4 sin(A0, )Af (£2, ) , 

(2.105) 

/„ = /„, + A -2A„,A r cos(A# )A/ (£2, ) , 

(2.106) 

A,, =/,„+/,. -2^,4 sin(A0,)M (£2, ) , 

(2.107) 

then by subtraction 


/„ - / 3 , = 4A ol A r cos(A0 f )M (£2,) 

(2.108) 

and 


/ 21 ~ / 4 , = sin(A0,)M (£2,) . 

(2.109) 


Using a similar method as was used for the static measurements and using the relation 


M(Q,W 0 (|a|) (2110) 

the data image for the time-average frames will become 

(l„ - 1 3 ,) 2 + (I 2 , - l 4.) 2 = 16 A*, 2 A, 2 J„ 2 (|ai). (2- 1 1 

Here the subscript t has been dropped from the amplitude variable for simplicity. The 
objective is to determine the fringe locus function which is the argument of the Bessel 
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function. The preferable method involves modulating the phase of eithei the object 01 the 
reference beams sinusoidally at the same phase and frequency as the test structure. This 
uses the fact that it is possible to shift J 0 fringes similar to the method in which phase 
modulation shifts consinusoidal fringes in conventional double exposure holographic 
interferometry [3, 5]. Adding a modulation step, a vibrating minor with a known optical 
phase, into the process makes it possible to solve tor the argument of J 0 - Mathematically, 
this modulation step is accomplished by adding a known variable into the argument of the 
Bessel function denoted as B. The variable B is commonly known as the phasor bias and 
takes the form of Eq 2. 1 1 2 

M(Q n B) = J n (p-B\). ( 2 - 112 ) 

Generally the intensity distribution can be expressed as Eq. 2.1 13 

im = I* + Im.Jo 2 [n«-B] < 2 * 113 > 

Where I a , represents the intensity or the energy of the local, average background of the 
scattered light and I mt is the local, maximum intensity of the scattered lighted. Equation 
2.1 1 1 a special case of Eq. 2.1 13 as demonstrated by the following equations 

In, = (I|t~ l3,) 2 + (12, - Ut) 2 ( 21 14) 

where the local, average background intensity is equal to zero in Eq. 2. 1 1 5 

I a , = 0, O' 115 ) 

and 

I mt =16A 0 2 A r 2 J„ 2 (|ai). ( 2 - 116 > 

The Bessel function is non-separable and in order to solve for Q t one should estimate the 
J 0 2 (|Q t |) intensity distribution as a cosine function, cos Qt [4].. This substitution can be 
made because the correlation between the crossings of the x — axis, the cosine function 
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and the Bessel function nearly coincide [3], Recall that a black fringe can be represented 
by every touching of the x-axis by the J 0 2 (|£2 t |) function. The general equation can now be 


rewritten as Eq. 2.117 


Ini fii I ml COS ' B], 

(2.1 17) 

Using the cosine half-angle theorem, cos~ v = 

-t(l + cos 2.v) Eq. 2.117 becomes Eq. 2.118 

In = la + Im-J ( 1 + COS[2fl - 2B] ) 

(2.1 18) 

Multiplying the cosine half angle through Eq. 

2.1 17 the new solution becomes Eq. 2.1 18 

In = l a + ^ Im + ^ Im COS[2f2 - 2B]), 

(2.119) 

where 


hi — hi ~~ hn, 

2 

(2.120) 

and 


Im ~~ hn. 
2 

(2.121) 

So, Eq. 2.118 reduces to Eq. 2.122 


Ih = I a + I m cos[2Q - 2B] 

(2.122) 

Where the bias modulation B is maintained at the same frequency as the object vibration. 

With three known values of B, (0, B and -B) Eq. 2.122 can be solved for the fringe locus 

function as follows: 


Ini = la + LC0S2Q, 

(2.123) 

I n 2 = l a + I m cos[2Q - 2B], 

(2.124) 

ho = I a + I m cos[2Q + 2B], 

(2.125) 
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Where l„i l „2 and I n j are die intensities of EOH images, each with a different phasor bias 
value. Using the cosine addition and subti action fonnulas, 

cos( x - v ) = cos( x ) cos( y ) + s i n ( v ) s i n ( r : ) (2.126) 

and 

cos(.v + y) = cos(.v)cos( v) - sin(.v)sin( v) (2. 127) 

the three values of B can effectively be used to solve for £2,. First, solving for the average 


background intensity I a from Eq. 2. 123 gives 

la = ini - Im COS 2£2. (2.128) 

Plugging Eq. 2.128 into Eq. 2.124 yields 

I n2 = 1„|- l m cos 2Q + I m [cos 2Q cos 2B + sin 2£2 sin 2B] (2. 129) 

Combining like tenns gives Eq. 2. 1 30 and then Eq. 2.131 

l n2 = l n) + l m [cos 2Q cos 2B + sin 20 sin 2B -cos 2Q] (2.130) 

then 

I n2 = I nl + l m [(cos 2B -1 )cos 20 + sin 20. sin 2B] (2.131) 

Solving for the local maximum intensity I m gives the following equations 

In2 - Inl = Im [(cos 2B - 1 )cos 2£2 + sin 2Q sin 2B] (2.132) 

I n2 - I nl = - I m [(1 - cos 2B)cos 20 - sin 20 sin 2B] (2.133) 


, ill inl (2.134) 

(l-cos2S)cos2Q-sin 2Qsin 2B 

Substituting Eq. 2. 134 in Eq. 2.1 25 yields Eq. 2. 1 34 

('.I -AJcos2fl 

" 3 (l - cos 25)cos 2Q - sin 2Q sin 2B 

(j n] - 1 n2 )cosQ cos 2fl-sin 2Qsin 2B (2.135) 

+ (l-cos2fi)cos2£2-sin2Qsin2fi 
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Expansion of Eq. 2. 1 35 gives Eq. 2.1 36 


I„,( 1 -cos 2B)cos 2 12 - I„, sin 212 sin 2B I n ,( 1 -cos 2B)cos 212 
+ I n isin 212 sin 2B = - I„i I cos 212 + l„: cos 212 
+ l„i cos 12 cos 2B l„i sin212 sin2B I,,: cos 212 cos 2B 
+ 1„2 sin 212 sin 2B (2.136) 

Solving Eq. 2. 1 36 for 12 yields Eq. 2. 1 37 


12 (a , y 


= — tan 

2 


(I- cos 28 )(/„., 

(sin 25X- 2/„, +/„ 2 +/„-,) 


(2.137) 


If the three intensities l n ,, I n2 , I„ 3 corresponding to three time-average EOH images with 
different phasor bias are substituted into Eq. 2.137 the result gives the approximate 
solution of the fringe locus function. Recall that the cosine function is merely an 
approximation used to solve for the argument of the Bessel function and though the two 
functions are similar they are not exact. Thus, Eq. 2.137 is only an approximate solution 


for 12, where 12, can be expressed as 

12, - 12, approx +e ( 2 ' 138) 

where e is the difference between the true Bessel function solution and the approximated 
cosine solution. Figure 2.8 is the plot ot ) solution versus the cosine solution, 

where the x-axis is the optical phase (fringe locus function) and the y-axis is the 
normalized fringe intensity [4]. Here it can be seen that the cosine function and the Bessel 
function do not match up exactly, which illustrates Eq. 2.137 is only the approximate 
solution. It should be noted that as the magnitude of the optical phase increases, the 
Bessel function and the cosine function more closely match up. It is for this reason that 
the squared Bessel function can be approximated by the cosine function squared. 
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Figure 2.8: Bessel function vs. Cosine Function. 

To correct for the error, £, a Bessel function look up table is created. The look up 
table is a table of the approximate cosine solution versus the exact Bessel function 
solution and is stored as a two dimensional array in computer memory [4]. The purpose 
of the table is to save computational time while the computer is processing the 
displacement data. Rather than spend the time computing the equations for each 
displacement condition, the computer will simply compute the approximate solution and 
then reference a preexisting table to find the exact solution. 

After computing fi ap prox, the holographic image may look something like Fig. 2.9. 
This image is known as a wrapped phase map, where the greyscale represent relative 
optical phase. The greyscale levels in this image have been scaled such that a greyscale 
level of 0 (black) represents an optical phase of (0 *tc), whereas a greyscale level of 255 
(white) represents an optical phase of (2*n). Thus the abrupt, discontinuous transitions 
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Figure 2.9: EOH wrapped phase map |4|. 

from black to white denote the occurrence of an EOH fringe, or additional multiple of 2k 
optical phase. Naturally the object surface does not have sharp discontinuities as implied 
by the wrapped phase map. The 0 to 2k (black to white) phase jumps are merely an 
artifact of the inverse tangent function used to compute the phase. To compute the final 
vibration amplitude map, the wrapped phase map must be "unwrapped to remove the 2k 
discontinuities. Phase unwrapping is often the most difficult part of EOH data 
processing, and there are numerous phase unwrapping algorithms in the literature ranging 
in scope of applicability and complexity. Conceptually, unwrapping is done by scanning 
a row of pixels across the image and adding up all of the 2 k changes [4]. The 2k changes 
will be positive or negative depending on whether the change on the wrapped phase map 
has gone from black to white or white to black. For example, if a single row of pixels is 
extracted from the center of Fig. 2.9 the line plot of Fig. 2.10 would be produced [4]. 
Here, each black peak represents a phase change of 2k. These changes are added or 
subtracted together to produce the unwrapped phase map (the red line). Finally, by 
dividing the unwrapped phase by the wave number, k — 2k/ A., and the sensitivity vector 
K of equation 2.37, the displacement is found. 





CHAPTER 3.0 


FRINGE CONTROL 


Electro-Optic Holography (EOH) is laser-based, diagnostic technique used to 
measure micron-sized displacements of static and/or dynamically loaded structures. 
Dynamically loaded structures are of greatest interest for this research project, where the 
displacement of the structure is great enough to cause the interferogram to become 
densely compacted with interference fringes. When the fringe density exceeds the ability 
of the software to clearly identify individual fringes the upper measuiement limit of the 
technique is compromised. By modulating the current source of the laser diode, the 
number of fringes can be reduced as a function of the modulation order. This also 
changes the intensity distribution as to that of higher order Bessel functions. Thus, by 
increasing the order, the sensitivity of the EOH system is decreased allowing for larger 
displacement measurements to be made. The reduction of fringes is theoretically 
unlimited but physically limited by the frequency bandwidth of the signal generator. This 
EOH fringe control method is known as Frequency Translated Electro-Optic Holography 
(FTEOH). 

This chapter focuses on the techniques and development for FTEOH. Section 3.2 
introduces FTEOH, where section 3.3 derives the theory behind the technique. Section 
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3.4 presents the methodology for estimating fringe density, which is followed by section 
3.5, which explains the FTEOH fundamentals. 

3.1 NOMENCLATURE 

A,, Amplitude of the object beam, the laser beam directed onto the test object. 

A n Amplitude of sideband components. 

A, Amplitude of the reference beam, the laser beam directed into the CCD 

camera. 

A v Amplitude of frequency shift. 

c Speed of light in a vacuum. 

CCD Charged Couple Device camera. 

d Temporary Variable, whose quantity is equal to 

[el v 2 + ALrAvl - 2e n vALA v„ cos(fl - ) ]' 

e General Amplitude of displacement. 

E 0 The time dependent displacement of a point on a vibrating object. 

e G Amplitude of the test object during cyclic vibration. 

EOH Electro-Optic Holography. 

E r The time dependent displacement of the PZT. 

e r Amplitude of the PZT during vibration. 

F Object frequency in Hertz. 

f Frequency in radians. 

F„ The object light field, the beam of light that illuminates the test object. 

F r The reference light field, the beam ol light that is directed into the CCD 

camera. 

Fu The light field reflected off of the vibrating structure. 
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Frequency Translated Electro-Optic Holography. 

A geometrical factor with g 0 = 2 for displacements perpendicular to the object 
surface and with normal incidence. 

A geometrical factor with g u = 2 lor displacements perpendicular to the object 
surface and with normal incidence. 

The light intensity of the laser beam. 

n th order Bessel function of the first kind. 

Sensitivity Vector, the vectorial difference between the propagation vectors. 
The wave number whose magnitude is given by 2k/X. 

Optical path length for the object leg. 

Optical path length for the reference leg. 

Optical path difference between the object and reference legs | L„ -L,\. 

Number of fringes masked by modulation. 

Characteristic function. 

Integer value used to increase the frequency of modulation. 

Index of refraction of the cladding material for the fiber optic cable. 

Index of refraction of the core material for the fiber optic cable. 
Proportionality constant. 

integer value used to increase the frequency of modulation. 

Remaining number of fringes after modulation. 

Exposure time of the CCD camera. 

Time in seconds. 

Total number of fringes without modulation. 

Temporary variable whose value is equal to e tt vs\n(ox + £2) . 
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Temporary variable whose value is equal to AiA ( sin( (iX +(ft v ) . 

Temporary variable whose value is equal to AvL, si n{ax +</>„). 

Frequency shift oflaser diode. 

Dirac delta function. 

Optical phase of the object beam. 

Optical phase of the reference beam. 

Unknown Phase, induced by modulation. 

Fringe multiplication factor. 

The wavelength of the laser light used for interference. 

The max angle beyond which rays that enter the fiber are no longer confined 
within the fiber. 

The critical angle beyond which rays will be bounded within the fibei. 

The approximate omega solution, made by approximating the Bessel function 
as a cosine function. 

The fringe locus function. 

Distance from origin to the maximum of the first peak of the J n curve. 

The step function for intensity jumps, which is equal to 271 radians. 

The exact fringe locust solution. 

The approximated omega with the addition of the step function. 

The unwrapped omega for the FTEOH technique. 

Angular frequency in radians of the test object. 


co Laser diode emission frequency, rad/s. 

v Laser diode emission frequency, Hz. 
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3.2 INTRODUCTION TO FTEOH 


The image of Fig. 3.1 illustrates a typical dynamic EOH fringe pattern, whereas 
the image of Fig. 3.2 illustrates how laser diode cunent modulation lias reduced the 
number of fringes. The objective ot this thesis is to develop a fringe conti ol method that 
will allow the user to represent the displacement field ot Fig. 3.1 with a fringe pattern like 
that shown in Fig. 3.2. 

Techniques to perform this type of fringe control have been developed for 
holography. For example, Aleksoff [5] presented a conventional film based holographic 
method that involved modulation of the reference beam using a mirror mounted on a PZT 
stack as means of controlling fringe characteristics. The mirror was forced at harmonics 
of the object vibration causing the frequency of the reference beam to be modulated. He 
stated that by modulating the reference wave any number of fringes corresponding to 
smaller displacements could be eliminated, thereby alleviating the holographic image of 
high fringe density. Although this type of fringe control has been applied to conventional 
holographic interferometry, fringe reduction has not yet been established for the EOH 
technique. 

The ideas of laser modulation presented by Aleksoff have been expanded upon in 
this thesis study and applied to the EOH technique. This new technique for obtaining 
large amplitude displacements has been given the name Frequency Translated Electro- 
Optic Holography (FTEOH). For the sake of definition it is best that the term fringe 
reduction be understood as it is applied to FTEOH. By masking the group ot fringes that 
represent areas on the test structure undergoing the small amplitude displacement (i.e. 
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those fringes nearest the zero order fringe) the image can be phase unwrapped and the 
displacement can be realized. However, the displacement that each fringe represents can 

A A 

no longer be approximated as — but now must be approximated as // ^ wneie /j is a 

fringe multiplication factor and is related to the order 01 the haimonic of frequency 
modulation of the laser diode current source. 



Figure 3.1 : Typical dynamic EOH image Figure 3.2: The desired result of 

of a cantilever beam. modulation. The density of the fringes 

from Fig. 3.1 have been reduced. 

Prior to the development of FTEOH many researchers [19, 28-30] used a 
technique in which the location of the zero order fringe was manipulated as means for 
determining large amplitude displacements. This method, known as phase mapping, is a 
relatively new technique which involves laser diode wavelength modulation as a means 
to determine the vibration amplitude of large displacements. Phase mapping heterodynes 
an external sinusoidal signal with the laser diode current source to produce a new current 

signal that in turn modulates the original laser beam and hence the phase of both the 

reference and object beams. Adjusting the amplitude and phase of an external signal 
controls the order of desired modulation. As a result, the location of the first bright fringe 
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is shifted from the nodal position to those points that vibrate at the same amplitude and 
phase as the reference beam [16, 19, 25-27]. By adjusting the phase and the amplitude of 
the external modulating signal the zero order fringe is shifted across the EOH image 
making it possible to map the vibrational amplitudes that may otheiwise be too great to 

be observed in the unmodulated case [25]. 

One of the big problems with the phase masking technique is that the zero order 
fringe must be driven to various positions on the object and more or less qualitatively 
reconstruct the vibrational amplitudes. Thus phase mapping, at least where image 
processing is concerned, is a non-automated process. This can be a significant pioblem 
in practical testing scenarios. Whereas the FTEOH technique has automated processing 
capabilities. 

From this phase mapping technique a fringe, FTEOH was developed. F I EOH 
differs from phase mapping by utilizing the frequency variation ot the external signal 
rather than the phase manipulation of the signal. The frequency of the signal is varied to 
produce an image with fewer fringes. These images have an intensity distribution related 
to higher order Bessel functions, which produces images where the fringes related to 
areas on the structure undergoing small displacement are masked as a result of 
modulation. 

The verification of the FTEOH theory allows the useful range of the EOH system 
to be extended well beyond former limitations, thereby making it possible to generate 
displacement data for amplitudes previously too great to be measured by the EOH system 
but too small to be measured through other conventional methods. Therefore, the theory 
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and practice of the FTEOH technique closes the gap between established structural 
displacement analysis and EOH. 

3J3 FTEOH THEORY 

The underlying theory of the FTEOH technique is presented in this section. The 
focus of the presentation is on the effects generated by modulating the current to the laser 
diode. In this thesis, the term "modulation” refers to modulating the current, which 
results in a corresponding modulation in the laser emission frequency. Physically, the 
laser diode current is modulated using a function generatoi connected to the laser diode 
current source. 

3.3.1 FTEOH REFERENCE AND OBJECT BEAMS 

In FTEOH the laser output is split into two beams, the reference beam and the 
object beam. The object beam illuminates the test structure. Light scattered fiom the 
vibrating test structure for the unmodulated case can be expressed as shown in Eq. 3.1 
F v {x,y) = 4,exp[i0 o (x,.y,O]- 

In Eq. 3.1, the light of the reflected beam is spatially varying depending greatly on where 
on the test structure the light is reflected. A <} is the amplitude of the reflected light and 
is the phase [16]. The phase of the object beam is given as Eq. 3.2 

27Cg„E„(x,y,t) G. 21 

<t> a (x,y,t) = 1 v ’ 

2nlX = k is called the wave number, g„ is a geometrical factor with g () = 2 for 
displacements perpendicular to the object surface with normal incidence! 16], E n is the 
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time dependent displacement of the point of interest on the vibrating test object and is 
expressed as Eq. 3.3 

I£ n {x,y,t) = e u sin(ft* + £2(0) ( 3 - 3 ) 

where e,is the amplitude of vibration of the object, a) is the angular frequency of the 
object, and £2(0 the phase. The object beam can now be expressed as Eq. 3.4 


F v {x,y) = A„ exp 


An 

i — e. 


X 


sin( <ur + £2(0) 


(3.4) 


The same evaluation can be given for the reference beam as shown by the following 
equations. Eq. 3.5 represents the unmodulated reference beam 
F r (x,y) = A r exp (i<p r {x,yj)) 


where 4, and <p, are the amplitude and phase of the reference light respectively. 
The reference phase is then expressed as Eq. 3.6 


A(x.,, ,).*5SpfL (3.6) 

where E, is the time dependent amplitude of laser diode current modulation and is given 
by Eq. 3. 7 

E r (x,y,t) = e r sin(<uf + <ft r (x,y,t). ( 3 - 7 ) 

The variables e r and <p r are the maximum amplitude and phase of the laser diode current 
modulation for the reference beam [25]. The reference beam can now be expressed as Eq. 


3.8 


F r = A r exp 


4 u 

i — e r sin (ax + <p r ) 

X 


(3.8) 
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where the argument (x, y, t) has been dropped tor simplification. Eqs. 3.4 and 3.8 can be 
given as the Fourier expansion (frequency spectrum) [5, 30] as shown by Fqs. 3.9 and 
3.10 


/•:</) = Z 


, 4/r 
J ( — e 
" X " 


exp (inLl)S(f - nF 


(3.9) 


W) = Z 


4 K 

V.( T ,» 


cxpiinfi, )t>( / - nF) 


(3.10) 


where ( j ) indicates that the light beam equations are given in tenns of the liequeney 


4 n 


domain [5], co = 2 ttF , F is the object frequency in hertz and -e„) is the n order 


Bessel function of the first kind with modulation depth (— e ,) . It should be noted that 

A 

the Bessel function appears in the Fourier expansion of the object and reference beam 
equations. 


3.3.2 THEORY OF MODULATION EFFECTS 

The major point of interest in laser diode current modulation, as far as FTEOH is 
concerned, is that both the reference and object beams are modulated, which forces Eqs. 
3.4 and 3.7 to be expressed as Eqs. 3.1 1 and 3.12 [25, 26]. 

. 4/r 


K = A, e x P 


i — {e„t»sin(rw + Q(r)) + AvL o sm(n ox + 0„)\ 
c 


(3.11) 


and 


F r 


4 jz 

A r ex p i — AvL f sin(fl*yf + <t> v ) 
c 


(3.12) 
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In Eqs. 3.11 


and 3.12 the wavelength of light is defined as 


— = — from Eqs. 3.9 and 

c X 


3.10. L t and L, are the optical path lengths for the object and relerence beams, 
respectively. The constant is the amplitude ol the liequcncy shift induced by 

modulation, and <j> l} is the unknown phase caused by reflection ott ol the vibrating test 

structure [25]. 

The intensity of the resulting modulated fringe pattern is determined by adding 
the reflected scatted light of Eq. 3.11 and the reference light of Eq. 3.12, and then 
multiplying by the complex conjugate as shown in Eq. 3.13, where * denotes complex 
conjugate. 

/ = (F, + F r )(F v + F r )’ 0-13) 


Expanding Eq. 3.13 the intensity becomes Eq. 3.14 


( 


l - \A„ exp 


— {e a vsm(ajt + Q) + AvL o s'm(nQX + 0„)) 
c 


+ 4 . exp ; 


^-AvL r sin(nax + <p v ) 
c 


, 

( 

j4, ex p 

— i 
v . 


+ A r exra - ;' 


An 

c 

An 


(e u vsin(ax + £2) + AvL a sm(nox + $ u )) 


— (AvL r s'm(ncot + </>„)) 


which can be rewritten as Eq. 3.15 


(3.14) 



An 

An 

+ A r 

An 1) 

/ = U, 

exp? — X 

exp? — / 

exp; — Z 


c 

c 


c JJ 


•< A. 


An 

exp- ; — / 

c 


.An _ 

exp; — Z 

c 


+ A. 


An _ 

exp - 1 — Z 

c 


(3.16) 


where 


X - <?„L»sin(ct# + £2) , 


(3.17) 
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Y = A vL it sin(//6!* + <p v ) 


(3.18) 


and 


(3.19) 


4 n 

■ 4/r V 

An 

exp/ — X 

exp/ — ) 

exp -/ — / 

c - 

c 

c 


Z = A vL, sin(/HMf +(pj . 

The expansion of Eq. 3.15 yields hq. 3.20 

/ = a; + A; + A, A, 

L ( JL L JL L J 

4 K . 4 k 4yr , ^ ta\ 

+ /!,/(, exp/ — A exp-/ — Y exp / — Z . p.^U) 

1 c JL c JL c 

Substituting Eqs. 3.17 - 3.19 into the third and fourth teinis ot Eq. 3.20 is given as Eq. 
3.21 

r . . 

-(e o Vsm(0X + Q) + A v(L, - L, )sin(/;trf + &„)) 


I = A] + A l + A a A r exp/ 


4 n 

c 


+ A u A r exp/ 


4 n 

c 


(e„t>sin(rrf + il) + At >(L, - L } )s\n( >uot + (p„)) 


(3.21) 


Since the CCD camera is a photon integrating device, the camera will only record the 
average light intensity during the camera exposure time. As a result, Eq. 3.21 must be 
integrated over time in order to obtain the correct intensity relation. 

J ( ' Idt = A] + A) + i F r J o r F> + jF; j' F v dt (3-22) 

The final term of Eq. 3.22 is of primary importance. It governs the intensity 

characteristics of the FTEOH interference pattern and is expressed as Eq. 3.23. 

k-f'v" 


_L rT 

T 


= — A {} A r exp/ j — [e„t>sin(&* + £2) + — L t )sin {ncot + <P„)] 


(3.23) 
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Let the path length difference between the object leg and the reference leg be given as 
A/ = | /. - L, | . The intensity for the image can be expressed as lx). 3.24 


I=\ M t \ 2 = 


exp/j — [e„usin(raf + Ll)+ Ai;AAsin(/j<tf + (p„ 


(3.24) 


where Mr - d„(f2,) is the characteristic function that denotes the interference of the two 
fields as a result of sinusoidal motion [3], Using the trigonmeteric identities of Eqs. 3.25 


and 3.26 


sin(x + y) = sin x cos y + cos x sin y 


(3.25) 


cos(x - v) = cos x cos v + sin x sin v 


(3.26) 


let 

d 2 = e]v~ + Au : A Is -e„l>AltALcos(i2-«^,J, (3.27) 

4/T . 

where d 2 is the square of exponential term of Eq. 3.24 (excluding as shown in 


equation 3.28 

d 2 = [e 0 t»sin(ryf + Q) + AuALsin(/J6X + <p u )\ . 

The characteristic equation can now be rewritten as Eq. 3.29 

4k n 


(3.28) 


M 7 


= — f exp((-/ifttf)exp/ — d 

T * C 


dt. 


(3.29) 


T - 0 " V c ) 

Using the Bessel function identity of Eq. 3.30 


exp i(d sin ox) = ^7 m (rf)exp/(rntuf) ( 3 - 3 °) 

m=-°° 

the characteristic equation can be written in terms ot the Bessel function as shown in Eq. 


3.31 
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(3.31) 


M r - ^ exp i{-nox)j,\ — d exp /(//»/ 6* ) . 

lit -— <« V ^ J 

The expansion of Eq. 3.31 yields Eq. 3.32 

M , = ^ J m I — [e 2 v 2 + A v 2 A Is - lev A vAL cos(fi - <p„ )f 1 - 


— f exp[/(/?7 - n)cot\dt . (3.32) 

T J) 

The final expression of Eq. 3.32 will be equal to zero for all eases except when m n . 

that is to say that the final term of Eq. 3.32 behaves as the Dirac delta function as shown 
in Eq. 3.33 


— f exp \i{m — n)(i)t\U —S(m — n). (o.33) 

T 

Thus, Eq. 3.33 only exists when the modulating frequency and the structural frequency 
of the test object match. The intensity distribution of Eq. 3.24 can now be expressed as 
Eq. 3.34 


1 C ) 


(3.34) 


Increasing the frequency of the modulating signal by harmonics of the test 


structure frequency changes the intensity distribution of the time average interference 
pattern. The intensity pattern is no longer governed by the square of the zero order 
Bessel function. Instead the pattern is governed by the ri h order Bessel function whose 
initial maximum occurs for increasing arguments of the Bessel function as shown in the 


Fig. 3.3. 
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Figure 3.3: The initial maximum of the 10 th order Bessel function occurs at a 
larger argument with respect to the initial maximum of the zero order Bessel 

3.3.3 Bessel Functions and the Intensity Distribution 

An EOH phase map is essentially a contour or mode shape map of the deflections 
incurred by the object during excitation. Thus each fringe represents a positive or 
negative out-of-plane displacement. The large, bright fringe about the periphery of the 
mode shape of Fig. 3.4 is known as the zero order fringe and is analogous to the term 
node or nodal point, which is a location of zero displacement on a vibrating object. When 
the laser diode current source is modulated at a sufficiently high order, the zero order 
fringe is eliminated from the image. This is also true of all sufficiently small amplitude 
fringes. 

Laser diode current modulation acts as a displacement filter that filters out or 
masks those fringes that represent the areas on the test structure experiencing low 


Tenth order vs. Zero order 
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amplitude displacements. In Fig. 

3.5, page 72, all three illustrations 
represent the same cantilever beam 
experiencing the same sinusoidal 
vibration. Figure 3.5 (a) is the 
interference pattern produced by 
conventional EOH techniques, i.e. 
no fringe masking. Figure 3.5 (b) 
is an exaggerated side view of the 
cantilever beam in oscillation, and Fig. 3.5 (c) is an image of the interference pattern 
produced by modulating the current to the laser diode. The dotted horizontal line that 
crosses through Fig. 3.5 locates the zero order fringe, the nodal position for the mode of 
oscillation. Notice that the number of fringes has been reduced from the fringe pattern on 
the right and the zero order fringe has gone black as a result of modulation. Regardless, 
the displacement field for both interference patterns is the same. 

Each successive increase in modulation order forces the initial maximum of the 
Bessel function to occur at larger and larger operands. Thus, as the order of the frequency 
modulation is increased it takes a greater object displacement to cause the initial bright 
fringe to occur. For those displacements that are too small no substantial phase change is 
seen and the image reports these areas as a dark null that increases with successive 
modulation order. The fringes that remain within the image are representative of areas 
experiencing of large amplitude displacement. The first bright, higher order fringe is not 
a nodal point as it is in the unmodulated case, but an area on the object with a 



Figure 3.4: No modulation with an excitation 
amplitude of 2 volts and a total of 7 fringes. 
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(A) Large Displacement 

(B) Small Displacement 
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I Front 

\ Modulated 

\ (C) 

^Number of fringes 
has been reduced 


Figure 3.5: Results of using the FTEOH method on a cantilever beam with sinusoidal excitation. 
The left image is of the unmodulated case while the right is the modulated case. Note the making of 
the low order frinoes for the modulated rase. 


displacement large enough to cause a significant phase change [25, 26], This comparison 
can be seen in Fig. 3.6 and Fig. 3.7 the ./ 2 (ft) vs. J 5 2 (Q) and 7 5 2 (Q)plots respectively. 

In short, as the modulating frequency order is increased, the intensity distribution 
changes as a function of the n 1 * 1 order Bessel function where each successive order masks 


Fifth order Bessel vs Zero order Bessel 



Fifth order Bessel Function 



Figure 3.6: J (( vs. J 5 


Figure 3.7: Fifth order Bessel function. 





the lower order fringes. Larger structural displacements can now be measured where 
before frequency modulation the fringe density would have been loo high to lesolve 
individual fringes. Now with the masking of lower order fringes each remaining fringe is 
clearer than before modulation, allowing phase unwrapping to be completed. By 
modulating the laser diode current source the number of fringes that will be unwrapped is 
fewer than the unmodulated case. Therefore, in order to obtain the correct structural 
displacement the distance that each fringe would represent would increase with 
increasing modulation order. 


71 


o 

GO 

03 


o 

UJ 


03 

(J 


s, 

O 

J3 

p 

C 

o 

CJ 




T3 

O 


O 

CJ) 


£ 


0 

1 r*; 

TJ 


(0 

Hi 


u. 

_l 

LU 

(/) 

C0 

Hi 

CQ 1 

HI 

01 

o 

2 

ri- 

co 

col 


as 


CO 03 O 

2 o -J 

*3 Cl) 1/1 
03 ZT 
OC P C3 

rn ■ — — 

o nn 

" o <? 

LL] o 

(-U H 

0 o P 

.If .1 

1 £ ^ 

e 

9> £ .2 
.z P o 
^ c 

CO ^ 33 

a af 

P C GJ 
O *C °2 

s 

■g o CQ 

^ tn U. 

j *| -8 

■s s ° 
§ s ? 

EC •£ J 

ill 



72 


Figure 3.8: EOH image and its corresponding FTEOH image. Here a 4th order modulation was used to mask two fringes from the FTEOH image. 
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Figure 3.9: Cantilever beam shown as first the EOH image and then with corresponding FTEOH images. It can be seen that by increasing the modulation 
order the interference fringes are in fact masked by a dark null. 




J23 
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Figure 3.10: More FTEOH images. The images here show a greater number of fringes, yet as the order of modulation is increased the fringes become 
more resolvable. 


J20 J24 J27 
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the modulation order is increased the interference 


3.3.5 LOOK-UP TABLES 


Equation 3.35 demonstrates how the approximate omega solution can be 
determined by approximating the Bessel function with a cosine function, derivation 
shown in Section 2.5. 

(l -cos2ff)(/„, - /„ : ) (3.35) 

(sin 2B)(2/„, 

If the three intensities l nI , I„ 2 , U corresponding to three time-average EOH images with 
different phasor bias are substituted into Eq. 3.35 the result is the approximate solution of 
the fringe locus function. 

The tan" 1 function can take values between K and -n and Eq. 3.35 will exist 
between /r/2 and -nil. Thus, the results of Eq. 3.35 will exhibit a discontinuity every 
time the phase increases or decreases by In rad. These discontinuities will be resolved 
in order to obtain a continuous phase map of the displacement. This process is obtained 
by adding or subtracting 2tt rad (O v ,„) every tune a discontinuity is encountered. 

Therefore, when the computer searches the wrapped phase map it searches for those 
discontinuities in the wrapped phase map that jump from black to white. If a 
discontinuity is found the step function, Q . sU . p , is added or subtracted to the approximate 


, , 1 . 

= 2 tan 


solution of Eq. 3.35 to obtain the unwrapped solution of Eq. 3.36. 


^unwrap ^upprox ~ ^ step 


(3.36) 


The cosine function is only an approximation of the Bessel function. Thus, Eq. 
3.35 is the approximation for the exact solution £i t where Q, is expressed as: 


£2| £2 un wrap £ 


(3.37) 
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where £ is the difference between the true Bessel function and the approximated cosine 
solution. To correct for the error, £, a Bessel function look up table is cieated. The look 
up table is a table of the approximate cosine solution versus the exact Bessel function 
solution and is stored as a two dimensional array in the computer memory [4]. The exact 
Bessel function solution is obtained by substituting the phasoi bias B (0. -B. B) into Eq. 
2.1 16 to obtain Eqs. 3.38 - 3.40. 

/. = \6A;A;Jl{Q, -0) 
l,, = 16 ^, 2 J 0 2 (Q ( + B ) 

/„ = \bAlA;Jl{Cl,-B) 
where /-, 1,, and /„ are the images captured with a zero, positive and negative phasor bias 
at each of the four phase steps of 0°, 90°, 180°, and 270°. Replacing the approximated 
values of Eq. 3.35 with those of Eqs. 3.37 - 3.40 the unwrapped solution becomes that of 
Eq. 3.41 

1 — cos(2i?) Jp '(Q, ~ 7?)--7 0 (Q, + B) (3.41) 

sin(2B) 2- 

The purpose of the table is to save computational time while the computer is 
processing the displacement data. Rather than spend the time computing the equations for 
each displacement condition, the computer will simply compute the approximate solution 
and then reference a preexisting table to find the exact solution. The resulting plot of a 
typical correction lookup table (B = ;r/3) is shown in Fig. 3.12, on page 81. 


G =-tan“ 


(3.38) 

(3.39) 

(3.40) 
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3.3.5. 1 LQQK-UP TABLE FOR FTEOH 

The Bessel function look up table for FTEOl I must account for the order of 
harmonic used to obtain the frequency modulated image, for example, if the fifth ordei 
harmonic was used to obtain a masked fringe image, the unwrapped omega, , of 

Eq. 3.41 would need to be a function of In general the unwrapped omega for the 

FTEOH technique is given as Eq. 3.42 

_ l _,[ l-cos(2ff) ./;(Q,-fl)-./;(£2,)+ff (342) 

"""' T “ 2 tdn sin(2/i) 2 • ./; (£2, ) - J; (Q, - «) - ■)] (S2, + B ) J 

The result of plotting the unwrapped omega for the modulated case, Eq. 3.42, versus the 

exact solution £2, , is shown in Fig. 3.13. 

The plot of Q versus £2, will not provide the correct solution because 
O is approximated using the cosine function, which does not contain any 

unwrap h'l * r 

modulated components. Thus, the plot of the newly calculated look up table is shifted 
down with respect to the original £2 look up table. A comparison of the two plots is 

shown in Fig. 3.14 where both the vertical and horizontal axes are in units of n. This 
downward shift is a result of FTEOH modulation and correlates directly to the area of 
masked fringes in a raw speckle image. As was discussed throughout chapter 3, the 
FTEOH image will contain a large dark null where lower order fringes have been masked 
as a result of modulation. Therefore, when the image is unwrapped the computer will see 
the masked area as having no intensity change, this corresponds directly to the zero 
amplitude area between the origin and approximately 6.4 £2, in Fig. 3.13, on page 81. 
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In order to correct the downshift of Fig. 3.13, the curve will need to be shitted up 
a distance, ft, , which is the distance from zero to the maximum of the first initial peak 

of the ./'(ft) curve to the origin, as shown in Fig. 3.15. This shilt will allow the curve of 
Fig. 3.13 to nearly imitate the lookup table of Fig. 3.12. By applying the ft, , shift to the 


unwrapped omega, Eq. 3.43 will be obtained. 

O — Q + £2 (3.43) 

Graphically the addition of the shift will nearly overlay the ft,,,,,,.,,,,, curve with that 

of the ft as shown in Fig. 3.16. This relation is again shown in Fig. 3.17 for 

Jl Jl f w and 4, where J; is the typical EOH ft look up table curve. Here it 

can be seen that the addition of ft, to ft will approximate the 4 (ft) look up 

table curve. A list of ft, *s are given in Table 3.1 . 

The solution for the FT EOH unwrapped displacement will not be the exact 
solution with respect to the EOH unwrapped displacement but will be a close 
approximation to the actual solution, hence the difference between the ft lookup 

table and the ft„ mra; , fT lookup table. 

The last and final step is to obtain the structural displacement for the test 
structure. Using the equations derived in Chapter 2 the displacement can be determined 
using Eq. 3.44. 

< 3 - 44 > 

K 
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Figure 3.1 2: The look up table is plotted by graphing the exact fringe locus function vs. the 
approximate fringe locus function for the zero order Bessel function case. The axes are shown 

in units of K radians. 



Figure 3.13: Exact fringe locus function vs. the approximate fringe locus function for the fifth 
order Bessel function case. The axes are shown in units of it radians. 








Figure 3.16: The exact fringe locus function vs. the approximate solution tor the zero order 
and fifth order (+ il ( , ) cases. The axes are shown in units ot it radians. 



Figure 3.17: The overlay of J 15 2 (C1), J.o 2 ^), J* l (0) and J„ z (£l) look up tables. For this figure 


figure J.l /. lilt uvnirtj UI « iu v — n ' ' * lm\ l l 

indicates all look up tables where FT varies in value from 0 to 15 m steps of 5. Note that the J 0 ( ) oo 
up table is an approximation for the higher order look up tables. The axes are shown in un.ts of * rads. 





Table 3.1: The n h order Bessel function and the location of the maximum oi the first initial peak. 



Locus of 


Locus of 

n 

maximum 

n 

maximum 

5 

6.42 

23 

25.32 

6 

7.50 

24 

26.36 

7 

8.58 

25 

27.38 

8 

9.65 

26 

28.41 

9 

10.72 

27 

29.45 

10 

1 1.77 

28 

30.48 

1 1 

12.825 

29 

31.51 

12 

13.88 

30 

32.54 

13 

19.43 

31 

33.57 

14 

15.98 

32 

34.59 

15 

16.96 

33 

35.62 

16 

18.06 

34 

36.64 

17 

19.09 

35 

37.67 

18 

20.15 

36 

38.70 

19 

21.18 

37 

39.71 

20 

22.215 

38 

40.74 

21 

23.25 

39 

41.77 

22 

24.28 

40 

. 

42.79 
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3.4 ESTIMATING THE FRINGE DENSITY 


This section details the method for estimating fringe using the I ThOH technique. 
This is not a phase unwrapping technique but rather a simple method toi quickly 
estimating the vibration amplitude. The displacement held lor the FT1.0H inteiteience 
pattern can be estimated using a simple linear relation that lies within the family ol 
Bessel functions and their relation to fringe masking. The number of fringes masked is a 
function of the Bessel function order, and as will be shown, not every modulation order 
will reduce the number of fringes from the FTEOH image. Modulating the laser dtode 
current source at any arbitrary harmonic will introduce a small error to the estimated 
displacement. Choosing a harmonic in which the first dark fringe ot the higher oidei 
Bessel function approximately aligns with a dark fringe in the unmodulated or the zero 
order Bessel function case can help to reduce this error. The alignment can be seen by 
plotting J] versus the n"' order Bessel function and choosing the function, which gives 
the best fringe alignment. The number of fringes that will be masked by the modulation 
order can be determined by counting the number of times to the left ol the alignment the 
J 2 (Q) touches the x-axis. Remember that touching the x-axis by a Bessel function 
represents a single dark fringe in the interference pattern. For example, in Fig. 3.18 the 
first dark fringe of the thirteenth order modulation ( J f, ) aligns with the 6 lh dark fringe of 
the unmodulated case (J 2 ). Here it can be seen that 5 complete fringes will be removed 
and the 6 lh fringe of the J] case aligns with the 1 st dark fringe of the case. Choosing 
the 13 th harmonic as the modulation order guarantees that no more, or no less than 5 
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fringes will be removed from the image. This particular alignment procedure can be 
found at every third harmonic starting at the 5 th harmonic (i.e., 5. 8, 10, 13, 16, 19 ). 


0 2 — 
0 18 - 
0 16 ^ 
0.14 - 
0 12 - 


Area of masked 




First white 


o.i - 
0.08 - 
0 06 - 
0.04 - 
0 02 - 
0 - 


First black 
fringe of Jo 




Figure 3.18: J, 3 2 vs. J () 2 . The 1st dark fringe of J 13 2 is the 5th dark fringe 
of Jo 2 . 


3.4.1 THE “EVERY THIRD ORDER” RULE 

To show that the “every third order” rule is merely an approximation, the zeros ot 
the n' h order Bessel function are listed in Table 3.2 along with the corresponding zeros of 
the zero order Bessel function. The alignment of the zeros is not exact but is relatively 
close. By using every third harmonic a simple linear relationship develops between the 
number of masked fringes and the modulation order as shown in Eq. 3.45 

* , 12 ( 3 . 45 ) 

MF = — n + — v ’ 

3 3 

where n is the n' h harmonic of the object’s frequency and MF is the number of masked 
fringes masked as a result of modulation. Figures 3.19-3.21 are plots ot the zero order 
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Bessel function versus higher order Bessel functions shown as an example of how Eq. 
3 45 in Fig. 3.15 is compared to ./ 0 , here it can be seen that the alignment of the 

first dark fringe of ./, 2 6 occurs at the 7 th dark fringe ol and as predicted by Eq. 3.45 six 
white to black fringes have been masked. 


Table 3.2: J n Zeros and the corresponding J 0 zeros. 


Jn 

First zero of J„ 
(Value of Fringe 
Locus) 

The corresponding 
zero of Jo 

Difference 

5 

8.77 

8.653 

0.117 


12.225 

11.791 

0.434 

10 

14.478 

14.930 

-0.452 

BHMHHHMI 

17.801 

18.071 

-0.270 

wrwmummmmmm 

21.085 

21.213 

-0.128 

19 

24.352 

24.352 

0.000 

22 

27.570 

27.493 

0.077 

25 

30.780 

30.635 

0.145 

28 

33.875 

33.776 

0.099 

| 'll 1 1 

37.156 

36.918 

0.238 

34 

40.330 

40.058 

0.272 

37 

43.487 

43.199 

0.288 

40 

46.645 

46.341 

0.304 
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Figure 3.19: J,, 2 vs. J„ 2 . The 1st dark fringe of J, 6 2 aligns with 
the 7th dark frintie of J«. 2 . 



Figure 3.20: J, 9 2 vs. J t > 2 . The 1st dark fringe of J t9 2 aligns 
with the 1st dark fringe of J 0 2 . 








If a non-optimum harmonic is used, like was done in Fig. 3.22, the last masked fringe is 
not completely covered. Instead, the 1"' dark fringe of is found to sit two-thirds of the 
way between two dark fringes of the ./,; case, liven though the interference pattern of the 

FT EOF! method will produce a fringe that is slightly masked, the hinge may be 
incorrectly counted as a whole fringe. If this is the case the estimated vibration amplitude 
will be erroneously too large. Therefore to eliminate this error it is best to choose one of 
the optimum harmonics and not have to correct for fractions of fringes when estimating 

the displacement. 

3,4,2 THE “EVERY THIRD ORDER” RULE AND ESTIMATING THE 
DISPLACEMENT 

Although showing how the “every third order” rule works on plots of the Bessel 
function, it is beneficial, to show how this relation actually corresponds to the pictorial 
representations of Fig. 3.18 - 3.22. In Fig. 3.19-3.27 the test object was excited at 925 
Hz. The images here show successive increase in modulation order starting with no 

modulation Fig. 3.23 and increasing 
to the tenth order in Fig. 3.31. The 
first order of modulation is not 
represented due to an increase in the 
sensitivity level [5, 31, 32], which 
results in an increase in the number of 
fringes. The second order modulation 
actually brings the number of fringes 

Figure 3.23: No modulation, 7 total fringes. 
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back to the original count by masking 
one fringe, so the lirst order in which 
fringes have been noticeably masked 
does not occur until the third order. 

Third order modulation has 
eliminated two fringes as shown in 
Fig. 3.24 when compared to the 
unmodulated image of Fig.3.23. 
These fringes were masked from the 
image by means of the third order of 
modulation, however, two fringes 
have also been removed from figures 
3.25 and 3.26, which use the forth 
and fifth modulation orders. 
Nevertheless, the image of Fig. 3.23 


Figure 3.25: 4 lh modulation order, 5 remaining 
fringes. 



should be used for estimating 
displacement according to the “every 
third order” rule, which predicts the 
alignment of fringes. According to 
the plot of T 5 vs. J 0 , the first dark 
fringe of J 5 aligns with the third dark 
fringe of ./„ , making the fifth 
modulation order the optimal choice 


Figure 3.26: 5 ,h modulation order, 5 remaining 
fringes 
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Figure 3.27: 6" 1 modulation order, 4 remaining 
fringes. 



Figure 3.28: 7" 1 modulation order, 4 remaining 
fringes. 



over the third and forth modulation 
orders. This entire synopsis is repeated 
for Figs. 3.27-3.31 where the optimal 
modulation order is shown in Fig. 3.28 
and 3.31, which are the eighth and 
tenth modulation orders respectively. 
This order will more closely yield the 
correct displacement. 

By using the “every third 
order” rule the structural displacement 
of the object can be easily calculated. 
To determine the estimated 
displacement, once the proper 
modulation order has been chosen, the 
operator will need to count the number 
of fringes in the FTEOH image 
denoted to as RF (Remaining Fringes). 
The RF number is then be added to the 
MF number of Eq. 3.42 to determine 
the total fringe count ( TF=MF + RF) 
of the unmodulated case. Next, a 
fringe multiplication factor r), is 
determined by dividing RF into TF 


Figure 3.29: 8 ,h modulation order, 3 remaining 
fringes. 
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Figure 3.30: 9 1 " modulation order, 3 remaining Fi Sure 3.31: I0 ,h modulation order, 3 remaining 

fringes. frin K es - 


number n = , which is then multiplied by — . For the modulated case a single 

RF 2 k 

black to white fringe has a displacement greater than one fringe in the unmodulated case. 
Thus, the multiplication factor rj is used to determine the linear displacement the 
structure for the masked fringe image. For example, if a 13 th order modulated image 
contained 4 fringes ( RF) the fringe multiplication factor is 


MF = —(13) + — = 5 
3 3 


MF + RF TF 9 

rj = = = — = 2.25 

RF RF 4 


each fringe now has a displacement value equal to 2.25 


2k 


( 3-46) 
(3.47) 
However, this 


methodology only works when optimum modulation orders are used, partial fringes are 
not accounted for in Eqs. 3.46 and 3.47. 
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3.4.3 THEORY SUMMARY 

Modulating the laser diode drive current results in a laser emission frequency 
modulation proportional to the input signal. Sawtooth variations of the external signal 
cause the phase relationship of the scattered object light to change and interact differently 
with respect to the unmodulated case. As a result, the optical intensity of the EOH 
images change such that each successive frequency variation ( na > ) causes the fringe 
intensity distribution to vary the EOH image as the n lh order Bessel function. Graphically 
speaking, the zero order fringe and all other sequentially lower order fringes from n - 0, 
1 n -l are essentially masked by the initial null of the n"' order Bessel function shown 

as an example in Fig. 3.32 using J 0 2 and f w . According to this plot, the 5 ,h dark fringe 
that will appear in the unmodulated, J 0 \ image will be the same as the first daik fringe 
in the modulated, Jf 0 image. The first bright fringe seen in the modulated image is not a 



First black 
fringe of V 

Figure 3.32: J i3 2 vs. J y 2 . The 1st dark fringe of Jj 3 2 « the 5th dark 
fringe of J 0 2 . 
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nodal point but an area on the object with large enough displacement to cause a 
significant phase change between the object and reference beams to create the n'" order 

bright fringe. This phenomenon is again illustrated in Fig. 3.33 and 3.34 with a V 0 \ 
unmodulated image and a modulated image. Here it can be seen that five fringes 
have been masked in Fig. 3.34, thereby verifying the higher order Bessel function 


intensity distribution theory. 



Figure 3.33: No modulation, 7 total fringes. Figure 3.34: 10 th modulation order, 3 remaining 

fringes. 
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3 5 FUNDAMENTALS OF FTEOH 


The goal of this section is to provide insight into how the FTEOH technique 
works based on signal analysis. This section does this by discussing the object, 
reference, and interference light beams in terms of their frequency content. 

Figure 3.35 is an enlargement of a portion of the FTEOH system. The figure 
shows a camera lens capturing a portion of the light scattered by the test structure and 
recombining it with the reference beam through a beam splitter. A CCD camera then 
captures the combined beam. Figure 3.35 illustrates the frequency content and physical 
signal processing that occurs in the FTEOH technique. A second order modulation was 
chosen merely to explain the theory of the FTEOH technique. 

It is shown that the light scattered from the object has been spatially modulated by 
integer harmonics of the vibrating test structure centered about the laser carrier 
frequency. The analysis uses the fact that a coherent light beam reflected by a vibrating 
object is temporally and spatially modulated in a way that depends on the type of 
vibration [11]. This is shown by the ao, a t , a 2 and a 3 spectral lines denoted in the object 
beam plot of Fig. 3.35. The relative amplitudes of these harmonic components are only 
estimated in the figure and will actually vary. The important point is that independent of 
the amplitude a component exists at hannonic frequencies. 

The reference beam is shown to have a clean second harmonic modulation from 
the carrier frequency, a single sideband suppressed carrier. This is obtained by using the 
sawtooth modulation, which is discussed section 3.5.2. Serrodynmg (sawtooth 
modulation) inherently suppresses all harmonics of the current modulation as shown by 
the reference beam. 
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The CCD camera integrates the intensity of the recombined object and reference 
beams to produce the interference beam signal shown in Fig. 3.35. The process of 
capturing the recombined light using a CCD camera is the same as performing a cross 
correlation of the two signals. A cross correlation is an operation between two different, 
yet continuous time energy signals. Upon combination ot the two signals a third, unique 
signal is produced. The correlation between two continuous-time energy signal, such as 
the reflected object beam and the reference beam, given as F v (x,y) and F r (x,y ) 
respectively, will have a cross correlation given by Eq. 3.48 [33] 

(r)= [j,,(t)F r (t-T)di. ( 3 - 48 > 

The term (/ — r) is the time convolution of the reference beam, however the same result 
can be obtained if the convolution of the reflected object beam were taken instead [33]. 

Figure 3.35 illustrates that the diode laser current was modulated at twice the test 
object frequency; this is seen by the large single spike at the J 2 position in the reference 
beam spectral profile. In the cross correlation process the J 2 component of the reference 
beam combines with the a 2 component of the object beam. Because of the cross 
correlation effect, the CCD camera sees the correlated portions of the signals as a 
constant. Since the second harmonic was the correlated part of the two beams the 
resulting interference pattern will have a fringe pattern governed by the square of the 
second order Bessel function. An additional feature of the CCD camera that can be 
applied is shuttering. This is an area that requires more research but the idea is that 

shuttering can be used to apply a sine filter (sine = , where T is the exposure 

period) about the zero frequency component, which in turn suppresses the sidebands of 
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the zero frequency signal content even further. It should be noted that in Figs. 3.8-3. 1 1 
and all of the FTEOH images, that the original zero order fringe is still present. It is 
expected that by reducing the shuttering speed, which in turn would decrease the width of 
the sine function, the presence of the zero order fringe in FTEOH images could be 
reduced. The mathematics for the shuttering effect of Fig. 3.35 is discussed below, where 
the time averaged intensity of the captured image over the period of exposure is given by 
Eq. 3.49 [5], 

I =l(F r \t)F(t)„dt (349) 

J 4) 


Camera Signal 



Reference Beam SSSC 


Figure 3.35: Frequency content and physical signal processing that occurs in the FTEOH technique 
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In Eq. 3.49 / is the time averaged combined intensity ot the object and reference beams, 
F*,(t) is the complex conjugate of the reference light beam, b, /t), is the reflected object 
beam. By taking the Fourier transform of Eq. 3.49, where the Fourier transform of the 
intensity is denoted by /. Eq. 3.50, which is a product of the sine function, is obtained. 

i = X(f) T s^c Tf (350) 

where 

xi.f)=K{f)*FXf) (151) 

x {f) is the cross correlation of the complex conjugate of the reference beam with that 
of the ob ject beam [5], 

T , sin/gf (3.52) 

nf 

This relation says that as T increases the width of the sine function decreases and filters 
out the frequency components that are not centered at zero. This effect is shown in Fig. 
3.35 where the range of the desired frequency components converge to /= 0 [5] . If the 
limits of the intensity are evaluated as the period T approaches infinity, as shown in Eq. 
3.53 

i = \imi r< XJT) = x(0)S(f) (3 ' 53) 

where S(f) is the Dirac delta function, then the function only exists at the desired 
freequency. With a long exposure time the value of the cross-correlation spectrum at the 
origin is the only value that becomes applicable [5]. This is why the FTEOH images are 
related to the higher order Bessel functions. 
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3.5.1 SCATTERING OF THE OBJECT LIGHT 

This section and section 3.5.2 discuss the scattering of the object beam and how to 
take advantage of the spectral components that result from reflection through use ot the 
sawtooth function. Modulation through use of the sawtooth or ramp function is known as 
serrodyning, which will produce a signal with a single sideband suppressed carrier 
(SSSC), which is instrumental in producing an FTEOH image with clean resolvable 

fringes. 

The reflection of the coherent object beam with that of the test structuie is 
temporally and spatially modulated in a manner that depends on the type of vibration [5, 
31, 32]. In this particular study, the vibration of the object is ot a sinusoidal manner. 
When the light is reflected off of the object’s surface, the modulated terms can be 
expressed as (exp(inQ)8(f - nF) ). That is, the reflected light is phase (exp(twQ)) as 
well as frequency modulated S(f-nF) and can be expressed in terms of higher order 

Bessel functions, as previously shown [5, 31]. 

The spectral composition of the light scattered by the vibrating object will contain 

the carrier frequency of the laser ( co c ) with side bands at integer multiples of the object 

frequency (Q) c ±nCO), where 0) is the object frequency and n is the n' h harmonic of the 

object frequency. The modulated spectrum is shown in Fig. 3.36 [5, 32]. 

The scattered object beam contains several different frequency components, 
although in comparison to the earner frequency (GHz), the components (Hz - MHz) 
might be considered inconsequential but in fact these components are of major interest 
for phase mapping methodology [5, 25, 26, 31, 32]. Using traditional modulation 
methods (i.e. a PZT actuated mirror), a conventional holographic image would be 
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Figure 3.36: Fourier Transform of the phase 
modulated spectrum of reflected light off of a 
sinusoidally moving object |34|. 

obtained by modulating the reference beam at integer multiples ot the object frequency 
(i.e., (O ±nco) [5], As a result, the synchronous or correlated component of the scattered 

object light and the modulated reference light is no longer the carrier signal a> ( , but is at 
O) c ±nco, which is located at the order side band [5, 31, 32] as shown in Fig. 3.36 
[34], The interaction of the scattered object wave with the modulated reference wave 
produces an image whose intensity is related to higher order Bessel functions. The order 
of the Bessel function is directly related to the n' h harmonic of the structure s excitation 
frequency, the frequency for which the reference beam is modulated. This method has 
been used by several authors [5, 31, 32] to increase as well as decrease the sensitivity of 
conventional holography. 
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3.5.2 SERRODYNING 


The spectral composition of the light reflected off a vibrating object tor the 
frequency modulated case is similar to that shown in Fig. 3.36 but the carrier frequency 
contains a time dependent component and is characterized by w =a) t ±Avsm pax as 
shown in Fig 3.37 [34], Each successive side band of Fig. 3.37, like those shown in Fig. 
3.36, contain a component of harmonic object frequency where the central component is 


J 5 (m) 



(co c +/-Avsin5(o) 

Figure 3.37: The modulated spectral composition of 
the scattered object wave 

the desired, modulated signal. The remaining components are the result of reflection off 
of a sinusoidally moving object. Unlike traditional modulation, these side bands do not 
contribute to the final EOH image. Instead, the synchronous component of the scattered 
object beam and the reference beam will always be the central frequency component, like 
the one shown in Fig. 3.33, regardless of the order of modulation. By modulating the 
laser diode current source both the object and reference beams have an equivalent 
spectral composition (©, = ± Ausin pox ). The only difference between the two beams 
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is that the object beam contains products of reflection and the reference beam does not. 
Therefore, when the light of the two beams recombine back at the CCD camera, the cross 
correlation will occur using the central frequency component of both the scattered object 
beam and the reference beam because these are the synchronous components of both 
signals. To insure that the intensity distribution of the interference pattern is governed by 
a desired higher order Bessel function a sawtooth waveform is used to suppress any 
sideband components. These sideband components are a result of laser diode current 
modulation via a sine function. 


3.5.2.1 SERRODYN1NG AT THE SOURCE 


Supposing that the external modulating signal is a sine function, instead of a 


sawtooth waveform, with a frequency equal to that 
of the test object’s excitation frequency [25], the 
optical frequency of the laser diode will oscillate in a 
sinusoidal manner about the lasing frequency with a 
frequency shift of Au , as shown in Fig. 3.38. 

For the particular case of FTEOH, an 
external periodic signal, generated by a function 
generator, is applied to the laser diode current 
source. However, the external modulating signal, 
given that it is a sinusoidal, causes the signal to 
oscillate about the lasing frequency at a frequency 
equal to that of the external modulating signal. The 
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Figure 3.38: Frequency shift of the 
laser diode. The solid line frequency 
spike represents the original lasing 
frequency. The dotted line is the 
frequency shift Av , which oscillates 
about the lasing frequency and is 
expressed as Ai>sin(cot), the 
modulating frequency. Where the 
modulating frequency = the test 
structure excitation. However, the 
mean frequency over the period of 
oscillation for the frequency is still 
the lasing frequency V . 
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optical frequency change can be expressed as shown in Eq. 3.54, where Av is the 
amplitude of the frequency shift and <p v is the phase of the source fiequency modulation 

[25, 26], 

A u = At>sin(tuf + 0„) (3-54) 

A small problem exists at the central frequency of both the scattered object beam 
and the reference beam. When a sine wave is used to modulate the frequency shift of the 
laser diode, the shift will contain a time dependent component. There is no guarantee that 
the components of the two beams will match up exactly. The scatteied object beam 
component may contain the positive carrier signal w =w + Avsm pax , while the 

reference beam may contain the negative carrier signal co c = co, -Avsm pax . 

If the Fourier transform of function shown in Fig. 3.38 is taken, the signal would 
contain several different spectral components within the frequency range 
oj -co +AvunpOJt to (o t = (O - Av sin p(Ot , which is not acceptable and highly 

noticeable within the viewing FTEOH images when the frequency of the modulating 
signal is increased. To correct for this, the sinusoidal waveform of the modulating signal 
is replaced by a sawtooth wavefonn. The sawtooth wavefonn will bias the signal to the 
positive side suppressing all undesirable sideband components. Modulating the laser 
diode current source by means of a sawtooth wave is known as serrodyning. The 
amplitudes of the residual sidebands after serrodyning can be expressed by Eq. 3.55 [39]. 

A „ =— [exp(- i(flco- nCo))dt 
Ik j 

= ! {exp[-/(Q-«)]-l} ( 3 - 55 ) 

2K(Q-n ) l 
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In Eq. 3.55 ftns the modulation frequency and A„ is the amplitude of the sidebands. The 
result, given that the frequency is modulated in integer steps, is a single sideband 
suppressed carrier (SSSC) with all residual sideband amplitudes suppressed to zero [39], 
In terms of its Fourier series expansion, the sawtooth function can be given as Eq. 
3.56, which can then be simplified as the summation of Eq. 3.57. 


v = — f sinox + — sin2rwf + -sin3rut + ... 
n l 2 3 


2A^ 

TC n= \ 


-sm ox 
n 


(3.56) 


(3.57) 


In the above equations A is the amplitude of the modulating signal. By substituting the 
Fourier senes of Eq. 3.57 into Eqs. 3.1 1 and 3.12, the object and reference beams can 
now be expressed as shown in Eqs. 3.58 and 3.59. 


An 

F v = 4 exp ' — 


e a vs>\r\{(Ot + Q) + AvL, ,^-sin(o* + (p v ) 


=i n 


(3.58) 


F = A r cxp\i— AuZ, r ^]-sin((<af + ^,)[ ( 3 - 59 ) 

r { c „=i n J 

The Fourier transform of Eqs. 3.58 and 3.59 are given in Eqs. 3.60 and Eqs. 3.61 




— (c„t> + A vL„) 
c 


exp (m(Q + (p v ))S(j - nco)r 




J n | —At »L r 

c 


exp (in <p, )S(f - n co) 


(3.61) 


(3.61) 


Thus, because of the presences of the Dirac delta function the above equations only exist 
at no), giving rise to the single sideband at no) with complete suppression of all spurious 
components. 
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3.5.3 WAVEFORMS. VERIFYING THE SAWTOOTH THEORY 


Based on the cross correlation analysis described above, a sine waveform cannot 
be used to perforin laser diode current modulation because its spectral composition will 
contain more than one component. After the first harmonic, using the sine wave, the 
images lose their contrast and definition and it is virtually impossible to clearly identify 
individual fringes. This same phenomenon is also true ot the triangular waveform. 
However, if a sine, square or triangular waveform were used to perform laser diode 
current modulation, a synchronous component is not guaranteed and as a result the 
fringes are unclear and difficult to distinguish, which is the case in Figs. 3.39 - 3.42. In 
Fig. 3.39 a sine waveform was used to modulate the laser diode current source at the 4 
harmonic of the test structure’s frequency. It can be seen that the fringes are difficult to 
define and do not resemble the fringe pattern of Fig. 3.39, an image created by a ramp 
waveform. In Fig. 3.41 a square wave was used as the modulating waveform, although 
the image is clear and all the fringes are resolvable, there is no fringe reduction. 
Therefore it can be concluded that the square waveform does not allow fringe reduction 
for the FTEOH method. Next, a triangular waveform was used to create the modulated 
image of Fig. 3.42. Unlike the case of the square wave, the number of fringes has been 
reduced with respect to Fig. 3.39; however, individual fringe resolution is poor, which 
makes displacement analysis extremely difficult. Based on the information obtained for 
these four waveforms it can be concluded that the best results for the case of FTEOH 
fringe reduction will occur when using a sawtooth waveform for modulation of the laser 
diode current source. 
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Figure 3.39: Forth order modulation obtained with Figure 3.40: Forth order modulation obtained by 

ramp waveform fringe resolution is clean and using a sine waveform fringe resolution has been 

distinguishable. greatly reduced. 



Figure 3.41: Forth order modulation obtained with Figure 3.42: Forth order modulation obtained with a 

square waveform. The image shows no signs of fringe triangular waveform. The fringe resolution is 
reduction. degraded. 
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CHAPTER 4.0 


EXPERIMENTAL SET UP AND PROCEDURE 

This chapter details the experimental setup used to obtain the data for FTEOH 
analysis. It includes a discussion of each optical component and how laser light 
propagates through the optical fiber based system. Although the system described here is 
fiber based, the FTEOH technique will work in exactly the same manner on a mirror 
based system. A fiber-based system was chosen tor ease of optical component 

movement. 

4.1 THE EXPERIMENTAL SET UP 

The experiments used to perform fringe reduction for the FTEOH system were 
done on the optical set up illustrated in Fig. 4.1. The ray trace of Fig. 4.2 was drawn to 
help visualize the laser beam path for Fig. 4. 1 . The optical system consisted of, in order 

of the journey of the laser beam: 

■ 100 mW laser diode, discussed in Appendix B. 

■ A telescope assembly consisting of: 

• A collimator, anamorphic prism pair and two converging lenses, all used 
to collimate, compress and shape the light beam. 
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■ A Faraday isolator to prevent back reflection into the laser diode, a beam splitter 
used to create two beams (the object and the reference beam) from the original 
laser beam. 

• Fiber optic couplers, used to direct the laser light into the fiber optic cable. 

■ Phase stepper and PZT mounted mirror, used to change the optical path length of 
object and reference leg respectively. 

■ CCD camera, used to capture the interference patterns. 

These optical pieces, which are discussed in detail in this section, were an integral 

part of producing FTEOH interference patterns for this thesis study. 



Figure 4.1: FTEOH Optical set up 
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Figure 4.2: FTEOH optical path. 


4.1.1 TELESCOPE ASSEMBLY 

As shown in Fig. 4.2, the light emitted from the laser diode is immediately 
directed through the telescope assembly shown in Fig. 4.3. This assembly is a 
combination of lenses and prisms that collimate and contract the diameter of the 
elliptically divergent beam of the laser diode. The telescope assembly was placed into the 
optical set up to more efficiently pass the laser light through the aperture of the Faraday 

Isolator shown as the right most element in Fig. 4.3. 

The first element in the telescope assembly is the collimator, illustrated in Fig. 4.4 
as a dotted cylinder. The light emitted from the laser diode is not a collimated beam with 
a circular beam profile, as is the case with all other gas and solid-state lasers. Instead, 
light emitted from the laser diode is highly divergent, with an elliptical beam profile that 
expands very rapidly from the laser aperture. Light emitted from the laser diode is not a 
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Figure 4.3: Telescope assembly and Faraday Isolator. 


collimated beam but a highly divergent, elliptical beam that expands very rapidly from 
the aperture of the laser diode. To correct the divergence of the beam, a collimator is 
used. The purpose of the collimator is to optically correct for the divergence using lenses, 
thus achieving a more collimated less diverging beam of light. The objective is to focus 
this point far enough away that the beam, in the area of interest, appears to be a long 
narrow column of light. Once the light is ‘’’collimated" it is directed through other optics 
with greater efficiency. The second element in the telescope assembly is the anamorphic 
prism pair shown as a cross-section in Fig. 4.5. This optic reshapes the elliptical beam 
into a nearly circular beam by magnifying along only one meridian [42], The light from 
the laser diode is characteristically elliptical, which causes a problem when attempting to 
pass light through a small circular aperture. Thus, an anamorphic prism is used to 
magnify the elliptical beam along one meridian creating a nearly circular beam. The 
transmitted beam is displaced laterally as a result of magnification as shown in Fig. 4.5. 
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Anamorphic prism 



Collimator 57mm 



Figure 4.4: Cut-away of the Telescope 
assembly. 


Figure 4.5: Cross-section of the 
Anamorphic prism pair [1|. 


The third and final element in the telescope assembly is a pair of double convex 
converging lenses. These lenses act together to contract the diameter of the circular 
beam. The beam expander was used to contract the laser beam diameter to obtain 
maximum light throughput through the Farady isolator aperture. 

The first lens in the telescope assembly of Fig. 4.4 is a 38 mm converging lens. 
Light passing through this lens will converge at a focal point 38 mm from the lateral axis 
and then will begin to diverge. If the second lens is placed at a distance of 19 mm from 
the focal point of the 38 mm lens the beam the now divergent beam of Fig. 4.6 will regain 
its collimated characteristics but will have a diameter one half the size of the original 


beam. 


38mm 19mm 

H M-H 



D/2 


Figure 4.6: Ray trace of the contracted beam. 
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4.1.2 THE FARADAY ISOLATOR 


After the light has passed through the telescope assembly it is directed thiough the 
Faraday isolator, a magneto-optic, used to prevent reflection from re-entering the cavity 
of the laser diode. A Faraday isolator is a unidirectional device based on the Faraday 
effect [43]. That is, the light will only pass through the device in the forward direction. 
The Faraday isolator is a magneto-optic device used to prevent external reflections from 
re-entering the laser diode cavity. Any light re-entering the cavity will disrupt the laser 
emission frequency stability. 

The isolator is comprised of polarizers, an optical crystal, and strong magnets. 
Combining the magnets and optical crystal produces a Faraday rotator, which can be used 
to rotate the polarization of the light passing through it [43]. This phenomenon, 
discovered by Michael Faraday, was found to occur when an isotropic dielectiic, such as 
a crystal, is placed in a relatively strong magnetic field. The presence of the magnetic 
field causes the dielectric crystal to become optically active, which in turn causes a beam 
of light, polarized the direction of the field, passing through the crystal, to become rotated 

in the direction of the magnetic field [36]. 

Light passing through in the forward direction, as shown in Fig. 4.7, of the 
isolator will initially pass through the input polarizer and become polarized along the 

mill fjsQ & 

,n p ut . . 

Polarizer > Polanzer 

Direction of Propagation 

Figure 4.7: The three major components of the Faraday Isolator. In this figure the light is traveling 
through the isolator in the forward direction |42]. 



Output 
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vertical plane. Upon passing through the Faraday rotator, the plane of polarization will be 
rotated 45°. The output polarizer, which has been aligned 45° relative to the input 

polarizer, will then let the light pass through unimpeded [43]. 

The isolator was designed to prevent back reflections from re-entering the cavity 
of the laser diode. Any reflection into this sensitive chamber will cause the diode to lose 
emission stability. In order to prevent this from occurring the isolator is placed in the 
direct path of the emitting beam. Any reflected light traveling in the opposite direction 
must pass through the output polarizer first and become polarized at an angle of 45°, the 
light then passes through the Faraday rotator and experiences and additional 45° rotation. 
The light is now polarized along the horizontal plane and will not be able to pass through 
the input polarizer, which only allows light polarized in the vertical plane to pass through 
unimpeded [43], The case for light traveling in the reverse direction is shown in Fig. 4.8. 


< > ■ Faraday 1 / * 

< — » I Rotator J 

Input 

Polarizer < 

Direction of Propagation 
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Figure 4.8: The polarization of back reflected light traveling in the reverse direction through the 
isolator. Light traveling in this direction will not be able to pass through the input polarizer and 
hence will be permitted from re-entering the cavity of the laser diode |39|. 


4.1.3 BEAM SPLITTER 

After passing through the isolator the beam interacts with a pane of glass that 
behaves as a beam splitter. The light both reflects and transmits through the glass creating 
two beams, which are at 90° angles to each other. The transmitted beam has been denoted 
as the reference and the reflected beam the object leg in Fig. 4.2. The reference beam 




contains most of the laser energy while the reference beam is typically only a small 
fraction, 5%, of the laser energy [4], 


4.1.4 FIBER OPTIC COUPLERS AND CABLE 

From the beam splitter the object and reference legs are directed along separate 
paths. Both legs, however, are eventually directed to fiber optic couplers. The tibei optic 
coupler has an optical lens placed at the head of a fiber optic cable. The fiber optic 
coupler receives the collimated laser light and focuses it down to a spot at the fiber optic 
core. The light successful focuses onto the core can then propagate through the fiber. 
Coupling the light into the fiber is a little more difficult than simply aligning the light 
through a small aperture. Figure 4.9 and 4.10 show both the front and the back of the 
fiber coupler used to direct the light of the reference leg into the fiber optic cable. 



Figure 4.9: Front of reference leg fiber 
coupler. 


Figure 4.10: Back of reference leg fiber 
coupler. 
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The optical fiber is made from silicon dioxide (SiCT) and dopant materials used to 
control the refractive index. Most fibers are constructed with an inner core with an index 
of refraction n„ slightly higher than the surrounding material n, called the cladding [44], 
The diameter of the fiber is approximately one optical wavelength wide and if not 
perfectly aligned through the center of the fiber the light will be dispersed in to the 
cladding material [44], For example, the light propagating through the fiber cross section 
of Fig. 4.11 is shown at two different angles. 

In the first case the light is at too great of an 
angle and enters the cladding material 
preventing propagation along the core. In the 
second case the light is shown at the 
maximum angle at which the light rays will 
be confined within the core. Thus, it is 
extremely important to guarantee alignment 



9 n : The max angle beyond which rays 
that enter the fiber are no longer 
confined within the Fiber 

0^ The critical angle beyond which 
rays will be bounded with in the 
fiber 


into the fiber is correct. Alignment using Figure 4.11: Propagation of light through a 

fiber optic cable, n c is the index of refraction 

. , .. . , j. . • for the cladding and n„ is the index of 

both the vertical and horizontal directions is refraction for the core |44F 

important, however, the fibers tend to be more sensitive to the pan and tilt adjustment. 
One will know when alignment has been properly achieved when the maximum amount 


of light has propagated to the end of the fiber. Alignment into the fiber for this 
experiment was made easier with the addition of the fiber collimator and the photo diode. 
The fiber collimator, just like the collimator used to collimate the beam of the laser diode, 
focus the light into the fiber. The collimator has an index of refraction that changes from 
the periphery radially to the center of the lens, which causes light to bend differently 
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along the radius of the lens. Light not in direct alignment with the fiber collimator will 
bend into the cladding material thus preventing light to travel the length of the fiber. 

The second element that aids in alignment is the photodiode. A photo diode 
greatly increases the ease with which the light beam can be coupled into the fiber by 
providing a signal that reacts when light has traveled the length of the fiber. The photo 
diode is placed into an open loop circuit, which is attached to the end of the optical fiber 
and an oscilloscope. The oscilloscope monitors the voltage change that occurs as the 
photo diode comes in contact with light. Once the light has been roughly adjusted into 
the index graded lens in the x and y position, the pan and tilt are adjusted until the 
reading on the oscilloscope indicates the maximum amount of light has been coupled 
through the fiber. The fiber ends for both the object and reference legs are then 
disconnected from the photodiode and placed back into the optical set-up. 

The light transgressing through the object leg travels through three meters of fiber 
before reaching the fiber end, as shown in Fig. 4.1 and 4.2. Once the light reaches the end 
of the fiber where it is no longer confined 
to the boundaries of the fiber core it begins 
to diverge very rapidly. The diverging light 
travels an additional meter through air upon 
reaching and illuminating the test structure 
of Fig. 4.12. This fiber end is placed at 
such a distance that the divergent light 
emitted from the fiber illuminates the entire 

Figure 4.12: Test Object. 

face of the object. This type of illumination 
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gives the operator the ability to choose a localized area or the entire surface ot the object 
to study. 

The structure, shown in Fig. 4.12, is an aluminum plate confined by a heavy 
frame and fixed at the base. The frame was used to increase the stiffness of the structure. 
On the backside of the structure a piezoelectric (PZT) actuator was placed near the 
centroid and was the source of dynamic excitation during testing. 

4.1.5 PHASE STEPPER 

The light traveling through the reference leg of the FTEOH system travels 
through approximately one meter ot fiber before encountering the phase stepper shown in 
Fig. 4.13 and 4.14. The phase stepper is composed ot a one inch diameter by one inch 
high piezoelectric cylinder, fiber optic cable and electrical leads with which voltage is 
supplied to the cylinder and cause the cylinder to expand. Approximately one meter of 
fiber was wrapped around the piezoelectric cylinder. Thus when the cylinder expands so 



Figure 4.14: Side view of the Phase stepper 


too does the fiber. The cylinder is excited by a stepped voltage, which causes the cylinder 
to expand and increase the optical path length of the reference leg. This expansion will 
travels and two additional meters before reaching the fiber end where the light is 
projected onto a beam splitter cube and directed onto the CCD array as shown in Fig 
4 . 15 . 



Figure 4.15: Reference light directed into the CCD 

4.1 .6 CCD CAMERA 

When the reference light is directed onto the beam cube splitter and recombined 
with the scattered object light speckle will be attained, ft speckle is not seen the 
interference fringes will also not been seen. Thus, it is vitally important that proper 
adjustment of the reference leg be attained. The speckle size of the speckles is also plays 
a major role in the unwrapping process. Specular interference (speckle) will result when 
the object and reference beams are combined at the beam splitter are the conditions of 
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mutual coherence and alignments are properly met. The observation ot speckle is a 
prerequisite of achieving EOH and/or FTEOH interferograms. 

Speckle size is also an important factor in the generation of interferograms. 
Ideally the speckle size should be approximately equal to or slightly larger than the 
physical size of a pixel on the CCD array. If the speckle size is too small, a single pixel 
may collect photons from multiple speckles whose interference may be out of phase. 
This results in an averaging effect that reduces interferogram contrast. 

Speckle size can be estimated by Eq. 4. 1 

DS = 2.44/1—, 

# 

where DS is the diameter of the speckle, A is the laser wavelength and — is the 1- 

number. The lens f-number is the ratio of the lens focal length divided by the lens 
aperture diameter [42], It is best to have speckles about the size of a single pixel, which 
can be controlled by adjusting the aperture of the camera lens. 

4.2 PROCEDURE 

Once proper alignment of the optical elements has been attained, acquisition of 
the data can begin. The lights in the laboratory were lowered to minimize background 
light collect by the CCD camera. The signal generator was turned on and the PZT, was 
excited at a fundamental frequency of the test structure. The amplitude of excitation was 
such that it induced a high fringe density. Next, the external signal connected to the laser 
diode current source, was excited at integer multiples of the structure’s frequency until 
the number of fringes reduced from the FTEOH image was satisfactory to the operator. 


With the phase stepper operating in four-step increments the data was captured and then 
unwrapped to yield the final results discussed in Chapter 5 Results and Conclusion. 


120 


CHAPTER 5.0 


RESULTS AND DISCUSSION 

In the previous chapters the theory of FTEOH was piesented. This chapter will 
detail the experiment and results that verify the FTEOH theory followed by a discussion 
on the future work for FTEOH. 


5.1 LIST OF EQUIPMENT 

Including the equipment discussed in chapter 4 the following components were 
used in obtaining the results for the FTEOH technique. 


• Infrared Laser diode: 

o Micro Laser Systems Lepton IV Series System 
o Model # L4852S- 1 1 5-TE-4 
o S/N 247 
o Power: 1 1 5 mW 
o Wavelength: 851 nm 
o Beam Diameter 

■ Horiz: 4.2 mm 

■ Vert.: 3.8 mm 


• Laser Diode Current Driver 

o ILX Lightwave Laser Diode Controller 
o Model # LCD-3724 
o S/N 37241228 
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• Function Generator 


o Hewlett Packer Multifunction Synthesizer 
o Model # 8904A DC-600 kHz 
o S/N 3518A08207 
o Opt8ZE, 002 


• Cantilever Beam 

o Length 7 in. 
o Width 2 in. 
o Depth 0.125 in. 

The cantilever beam used for analysis is shown in Fig. 5.1. The shaded area of Fig. 5.2 


represents the actual area of study. 



Figure 5.1 : Cantilever beam used to 
obtain EOH and FTEOH data analysis. 



Figure 5.2: Shaded area represents the 
area used to during data acquision of 
this study. 
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5.2 THE EXPERIMENT AND DATA CAPTURE 


Experimental verification of the FTEOH theory was accomplished using a 
vertically mounted cantilever beam. The beam was sinusoidally excited using a PZT 
mounted at the beam center near the tree end. The cantilever beam was excited at 428 Hz 
with amplitude of 3.30 V. In this test setup only a portion ol the cantilever beam of Fig. 

5.1 was analyzed as illustrated in Fig. 5.2. 

Data was captured in two sets. The first set was acquired using conventional 
EOH techniques, i.e. no laser diode current modulation. The second set was acquired 
using FTEOH techniques with sawtooth modulation applied at 2140 Hz, the fifth 
harmonic of the test structure excitation. The test structure excitation amplitude was 
constant for both cases. 

In previous test cases the phasor bias settings were determined separately for the 
modulated case. The process of determining these settings is not as easily followed as in 
conventional EOH, mainly because the peak of the first fringe of a higher order Bessel 
function is more difficult to identify when the object excitation is turned off. This is a 
result of a lower peak intensity for the ever increasing Bessel function. It was also 
determined in those cases where the phasor bias settings were calculated precisely, that 
they were nearly the same as the setting for the conventional EOH system. Thus, it was 
decided to approximate the amplitude for the phasor bias to that of the conventional EOH 
phasor bias amplitude the results would produce the same displacement. 

Here it can be seen that two fringes have been masked from the FTEOH image 
leaving behind a dark null. Although less apparent, the two masked fringes can be seen 
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in Fig. 5.4, these effects were discussed in Section 3.3. The procedure used for capturing 
data for both the EOH and FTEOH technique is detailed in Table 5.3. Figures 5.3 and 
5.4 show the interferograms viewed for conventional EOH and FTEOH , lespectively. 


Table 5.1: Data capture procedure for both the EOH and FTEOH techniques. 


step 

Procedure for data capture 

i 

The First step is to determine the amplitude (voltage) at which the phase stepper 
will be operating to produce the 0°, 90°, 180° and 270° phase shift required in Eqs. 
3.35-3.37. 

2 

Produce a static interference pattern. 

3 

Slowly increase the amplitude ot the phase stepper, from zero, until the static fringe 
pattern shifts one complete black to white and back to black fringe. Record this 
voltage. 

4 

Excite the test structure at the desired frequency and amplitude. Record the 
frequency, &), and amplitude. At this point several cleaily identifiable fringes 
should be present on the computer monitor 

5 

For the leg that contains the phasor bias, the PZT mounted mirror will need to be 
excited at, 0 ) , the frequency of the test structure. This can be accomplished by 
connecting the phasor bias (the PZT mounted mirror) to a function generator and 
dialing in the desired frequency. Set the amplitude at zero and slowly increase the 
amplitude of the phasor bias PZT, via the function generator, until the zero order 

fringe is lost. 

6 

Leaving the amplitude set, adjust the phase of the phasor bias excitation, via the 
function generator, until the zero order fringe is regained and its width is 
maximum. Record the phase, a . Be careful not to increase the amplitude in step 5 
past the point at which the zero order fringe is lost. If this point is significantly 
passed you will not be able to regain the zero order fringe by adjusting the phase. 

7 

Turn off the test structure’s excitation (0 Hz) and reduce the amplitude of the 
phasor bias (PZT mounted mirror) to zero. DO NOT CHANGE THE PHASE. The 
structure should appear completely white on the computer monitor. 
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8 

Slowly increase the amplitude of the phasor bias (via the function generator) from 
zero until the entire structure goes black. Record the amplitude, V/ . At this point 
the value of B, from Eq. 2.1 10, is equal to 2.405, the argument of the first dark 
fringe of Jo, as shown in Fig. 5.5. 

9 

The voltage K/ v ' will need to be reset by a scaling factor. This factor is the desired 
magnitude of the bias vibration B divided by 2.405. The scaling factor is denoted 
as for bias voltage and can be computed by using Eq. 5.1 

V = B V (5-1) 

B 2.4045 

In Eq. 5.1 the desired value of B for this thesis was ^ , however any known value 

can be used, as indicated by Eqs. 2.121 - 2.123. In Chapter 2 it was shown that 
Eq. 2.120 contained three unknowns, /„. and 12, thus, by adding three known 

values of B into Eq. 2.120 there will now be three equations and three unknowns, a 
solvable system of equations. Following this reasoning three sets of images with 
three different, known phasor bias settings (0, B, and B) are captured at each of 
the four phase angles (0°, 90°, 1 80% and 270° of Eqs. 2.98 - 2. 1 0 1 ). Therefore, Eq. 
2.120 can now be solved for the fringe locus as shown in Eqs. 2.120 - 2.133. 
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1 


Square of the zero order Bessel function 



Figure 5.5: Zero order Bessel function. 


5.3 FTEOH Data Capture Procedure 

When taking FTEOH measurements, the amplitude of the laser diode current 
modulation must first be calibrated to ensure clear resolution of the individual fringes. 
This procedure is outlined in Table 5.2. If the amplitude is not set correctly fringe 


masking will not occur. 

Table 5.2: Setting the amplitude of the laser diode current source. 


Step 

Setting the Amplitude of the Laser Diode Current Modulation 

1 

Tune the laser diode current modulation and the test structure to the same 
frequency. 

2 

Turn the excitation of the test structure off. 
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Starting with the amplitude of the current modulation at zero, slowly increase the 
amplitude until the image of the test structure goes black. This is the calibrated 
amplitude for obtaining FTEOH measurements and should be used exclusively 
throughout the measurement cycle. At the point where the test structure goes 
black is the first null of the zero order Bessel function. 


5.4 DETERMINING DISPLACEMENTS 

After data capture the EOH and FTEOH images were phase unwrapped in 
accordance with the procedure of section 3.3.5. As a part of the unwrapping process the 
image (EOH and/or FTEOH) is dissected into small columns, where omega is then 
determined for each column. The three-dimensional plot for both the non-modulated and 
the modulated case are shown in Fig. 5.6 and Fig. 5.7. The rows of the columns of Fig. 
5.6 and Fig. 5.7 are then averaged across the plot. The average value of each row is then 
divided by the sensitivity vector, Eq. 3.41. The overall result is the plot of average row 
value versus the displacement of the test structure. Figure 5.8 represents the average row 
value versus the displacement for both the modulated case and the non-modulated case. 
Here it can be seen that for the point of interest, which is the base of Fig. 5.3 and Fig. 5.4, 
on the test structure the overall displacement for both cases will result in the same value. 

In Fig. 5.8 the displacement immediately after the dark null does not match with 
that of the non-modulated case this is most likely a result of the non-linearities of the first 
bright fringe of the modulated case. The square of the zero order Bessel function can be 
approximated as a cosine wave where the approximation becomes more exact as the 
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Displacement {m\ 



Figure 5.6: 3-D Plot of displacement for the interference pattern of Figure 5.3, i.e. no modulation. 
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Displacement im) 



Figure 5.7: 3-D Plot of displacement beam and interference pattern of Figure 5.4. 
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operand of the function increases. The same is true for higher order Bessel functions, the 
cosinusiodal approximation between the two functions becomes better and better as the 
operand increases. Initially the approximation does not match well because the width of 
the first couple waves is wider than that of a cosine wave ol the same frequency and thus 
the cosine approximation used when determining £2, h as a greater error at for those 

fringes immediately after the null. 

Here it can be seen that the overall displacement for the test structure for both the 
modulated and the non-modulated is case is 5.7x10 7 m. This displacement can be 
verified using the displacement equation given by Vest [22], shown in Eq. 5.1 

n p) = fa • (5A) 

An! 

In Eq. 5.1 L(P) is the displacement of a particular point on the test structure, is the n 
zero of J 7 (Q), shown in Table 5.3. One could calculated the displacement for the 
interference pattern of Fig. 5.1 by counting the number of fringes ~5.5, determining the 
approximate operand just past the 5 Ih zero and using those values in Eq. 5. 1 . 


Table 5.3: Zeros of the Bessel function J« j 


n 

L 

n 


1 

2.2048 

ii 

33.7758 

2 

5.5200 

12 ~l 

36.9170 

3 

8.6537 

13 

40.0584 

4 

11.7915 

14 

43.1997 

5 


15 

46.3411 

6 


16 

49.4826 

7 

21.2116 

17 

52.6240 

8 

24.3524 

18 

55.7655 

9 

27.4934 

19 

58.9069 

10 

30.6346 

20 

62.0484 
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The operand for 5.5 fringes can be determined by averaging the distance between the 5 
and the 6 lh zero ~ 16.50095. Now using these values the displacement can be calculated 
using Eq. 5. 1 as shown below: 


L(P) = 


16.50095A 

8/r 


5.587x1 O' 7 /;/ 


(5.2) 


where A = 851x10'" m . The calculated value of Eq. 5.2 is very close to that of Fig. 5.7, in 
which case one can conclude that the displacement tor Fig. 5.1 and Fig. 5.2 is correct. 


5 5 FUTURE WORK 

Future work for this project would include the investigation of shutter control. 
Although discussed in this thesis, shuttering of the CCD camera was not applied to any of 
the data gathered. It is the belief of this author that if shuttering were applied to FTEOH 
the appearance of the residual zero order fringe would either be completely masked 
having no difference in intensity from that of the other masked fringes or have an 
intensity that is significantly darker than those images shown here. 

5 6 CONCLUSION 

The primary objective of this study was to show that modulating the frequency of 
the laser diode current source produced an interference image with masked lower order 
fringes and the displacement of the modulated image was the same as the non-modulated 
image. It has been shown that this is indeed the case and although it can be difficult in 
practice the FTEOH technique can be used to expand the useful measurement range of 
electro-optic holography. 
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